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SUMMARY 


The objective was to determine if a pressure release boundary condition can 
be achieved by coating an elastic shell with a visco-elastic material. One necessary 
condition is that the coating must acoustically decouple the shell from the 
scattering problem. Two closed cell rubbers and two cork-rubber composites 
(nitrile and neoprene based) were investigated. The dynamic viscoelastic constants 
of the materials were determined by wave propagation techniques. The far field 
scattering form functions for an infinite cylindrical shell coated with the viscoelastic 
material were calculated using the complete elastic equations of motion. The form 
functions were experimentally measured for the different materials at different 
thicknesses as verification of the theory. A thick finite right cylindrical shell was 
coated with .25 inches of closed cell neoprene and the normalized scattered 
pressure measured. The pressure release normalized scattered pressure was 
determined for the end on incident plane wave case using the acoustic radiation 
Simplified Helmholtz Integral Program (SHIP). The pressure release normalized 
scattered pressure was determined for the side incident case using a modified 


Combined Helmholtz Integral Equation Formulation (CHIEF) radiation program. 





The material property measurements showed the closed cell rubbers have 
longitudinal wave propagation speeds of approximately 150 m/sec and 
attenuations of 30 dB/cm. The cork-rubber composites have longitudinal wave 
speeds of approximately 300 m/sec and attenuations of 7 dB/cm. 

The scattering measurements demonstrated that a thin shell (inner radius to 
Outer radius ratio of .97) could be made to scatter in a pressure release manner 
with a .25 inches of nitrile. The rubber-cork composites could not produce the 
pressure release effect for nondimensionalized wave number (product of the wave 
number and the radius of the cylinder) values less than 4 with reasonable 
thicknesses. The coated finite thick shell, with side incidence, scattered within 1 
dB of the pressure release theory and completely eliminated the first 
antisymmetric Lamb mode circumferential wave. For the end on case, the coating 
significantly reduced the resonance effects and eliminated the scattered signal from 


the back end of the shell. 





CHAPTER I 


Introduction 


Acoustic scattering from objects has been studied extensively since the days 
of Lord Rayleigh. The most recent advances, which involve extensions of the basic 
theory to more arbitrary shapes and mixed boundary conditions, can be directly 
attributed to improved numerical techniques and the digital computer. There are 
two fundamental, degenerate boundary conditions which may occur on the surface 
of the scatterer. The first is the rigid surface or Neumann boundary condition 
where the normal acoustic velocity of the surface is zero (i.e. the density and sound 
speed of the scatterer + ), The second is the pressure release surface or Dirichlet 
boundary condition where the acoustic pressure is zero on the surface (i.e. the 
density and sound speed of the scatterer + 0 ). Reductions to these two degenerate 
conditions are traditionally made when the specific acoustic impedance of the 
scatterer is significantly greater than (for the rigid case) or less than (for the 
pressure release case) the specific acoustic impedance of the surrounding medium. 

The rigid body case (Neumann condition) and the pressure release case 
(Dirichlet condition) have been extensively researched analytically for simple 


shapes!. The rigid body has also been the subject of extensive experimental inves- 





tigation. Since a perfectly rigid body does not really exist, this condition may be 
considered a limiting form of the more precise problem of scattering from elastic 
objects where the degree of elasticity is specified or to be determined. Elastic body 
scattering is a current major research area and the review papers by Uberall2,3and 
Nigul4 contain nearly one hundred references on this subject. 

The pressure release surface (Dirichlet condition) has not received the 

eiiental attention of the rigid body case for two reasons: 

1. This condition is not observed in air. 

2. Until recently materials such as closed-cell rubbers and corprene which 
approximate a pressure release surface in water have not been 
commercially available. 

The use of the closed-cell rubber as a decoupling material has provided a practical 
application. Since the pressure release material is not structurally suitable as the 
scattering body itself, it will be considered as a surface coating. The rigorous 
theoretical approach to the coated scatterer is to apply the exact eases 
equations to the structural material (with or without absorption) and to apply the 
exact viscoelastic equations to the coating. The scattering problem is then solved 
as a boundary value problem. If the coating were truly a pressure release material, 
the elastic behavior of the structure could be neglected. This decoupling of the un- 
derlying elastic structure from the fluid allows the scattered field to be determined 
by treating the body as an impedance surface (using only the boundary conditions 


at the exterior surface) and solving the exterior field problem using well 





documented techniques. 

This thesis investigates the pressure release surface condition in detail. The 
procedure to be followed is that suggested by Neubauer. The backscattered 
pressure field for a pressure release, finite cylinder is determined using the 
modified Naval Ocean Systems Center (NOSC) CHIEF program for side incidence 
(or normal incidence) at low frequency and the Naval Research Laboratory (NRL) 
SHIP program for end on (axial) incidence. The backscattered pressure of a 
coated finite cylinder is measured and compared with the theoretical pressure 
release results. The theoretical prediction of the backscattered pressure for an 
infinite coated cylinder using the experimentally determined material constants is 
used for comparison with measurements on experimental models. Finite cylinder 
data has compared well with infinite cylinder calculations when long cylinder are 
used and the ends are either outside the main lobe of a piston type projector or the 
length of the cylinder is greater than the length of the first Fresnel zone, i.e. 
farfield of the cylinder, using a point source®, The overall objective is to determine 
to what extent the pressure release condition can be approached, (i.e. for what fre- 
quency range and for what thickness of coating material), with some commercially 
available materials. Finally, the significant question of when the scattering field 


can be determined treating the surface as a simple pressure release boundary or 





the more complex elastic-viscoelastic solution must be use in order to achieve the 


required accuracy is answered. 


Background Review 


Scattering From the Infinite Cylinder and Cylindrical Shell 


The first solutions to plane wave scattering by an infinite cylinder were 
developed by Lord Rayleigh’. He considered both rigid, immovable and nonrigid 
cylinders with diameters small compared to the acoustic wavelength in the fluid 
(ka< <1 where k = w/c, a is the radius of the cylinder, w is the angular frequency, 
and c is the speed of sound in the fluid). Rayleigh also outlined a method of 
solution for larger diameter cylinders and spheres using an infinite series of cylin- 
drical and spherical harmonics, respectively, which has become known as the 
Rayleigh series. The solution method using the Rayleigh series is called the 
normal mode solution method. Morse® demonstrates the normal mode solution 
method for ka values up to 2.5. In 1951 Faran? published a classical paper on the 
exact solutions to scattering of plane waves from an infinite, homogeneous, isotrop- 
ic cylinder capable of supporting both shear and compressional waves. His method 
of solution was to solve the elastic equations of motion developed by Love!9. The 


three dimensional equation of motion in vector form is 





-xcellent experimental confirmation for elastic spheres and spherical shells 11-16 


: ° ° 1] eye Sas 
dickling!1 was the first to utilize the digital computer for ka values less than 20 





and found convergence of the harmonic series to be extremely slow as ka 
increased. He also formulated the result in terms of a far field form function, f,, 
which normalizes the scattered pressure, P.(@), to the incident pressure, P,, radius 
of the scatterer, a, and the distance to the field point from the acoustical center of 
the scatterer, R in such a manner so as to be equal to 1.0 for the rigid body case. 


The cylindrical far field form function is 








(1-2) 


Neubauer et al.!5 demonstrated the first quantitative agreement between normal 
mode theory and experimental results for solid elastic spheres. Dragonette et al. 16 
demonstrated that good quantitative results could be obtained using short 
broadband pulses and analyzing the frequency content of the scattered signal. 
Experimental verification of the normal mode solution method for the infinite 
cylinder (using a finite cylinder) was obtained by Dardy et al.17, 

Due to differences in the scattered fields, it is necessary to differentiate 
between scattering from cylindrical shells and scattering from solid cylinders. 
Barnard! observed numerous periodic echoes for an incident pulse on solid and 
hollow brass cylinder which suggests additional wave contributions besides the 
purely specular reflection. The hollow cylinder (shell) was also observed to scatter 


differently than the solid cylinder. Subsequent authors!9,20,21 proposed the exis- 





tence of two types of circumferential waves for the solid cylinder. The two types of 
waves are: 1. Franz waves which are analogous to the geometrically diffracted cir- 
cumferential waves discovered by Franz22 while investigating diffraction of elec- 
tromagnetic waves by perfectly conducting spheres and cylinders. 2. Rayleigh type 
(or R type) waves which are surface waves which depend on the elastic properties 
of the scatterer. Neubauer? gives an excellent review of the literature and theory 
of circumferential waves which is briefly summarized below. 

The two wave types generated by the solid cylinder resulted from the 
observation that one wave, the Franz wave or geometrical wave, had a wave speed 
below the propagation speed of the fluid and the other wave, Rayleigh type waves, 
had a wave speed significantly above the sound speed of the fluid medium. 
Neubauer24 and Harbolt and Steinberg2> experimentally verified the presence of 
the geometrical or Franz waves. Theoretical evidence for the Rayleigh type waves 
was developed by Grace and Goodman? with experimental verification by 
Bunney et al.27 and Neubauer28, Figure 1-1 illustrates the the specular reflection 
(backscattered) contacting the cylinder at point A and returning. Also illustrated 
are two Franz waves originating at points B and C and traveling in the fluid around 
the back of the cylinder. The Rayleigh type waves travel paths similar to the Franz 
waves but Beict in the solid cylinder not in the fluid. 

Although the normal mode solution method leads to the correct scattered 
pressure verses frequency result, the series converges very slowly and becomes 


computationally impractical for ka values greater than 20. Another drawback of 
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this method is that it does not easily isolate the individual mechanisms, i.e. the 
specular and circumferential waves, contributing to the scattered pressure field. 
Uberall et al.2? has demonstrated the relationship between circumferential or 
creeping waves and the normal modes, but unfortunately the relationship is not 
Straight forward. The current method to extract the individual wave type 
contributions is the Sommerfeld-Watson transformation of the normal mode 
Series. This transform maps the normal mode series into the complex plane in the 
form of a contour integral. Two groups of poles are observed which yield the two 
types of circumferential waves. The Sommerfeld-Watson transformation has been 
extensively used. It converges rapidly and its results agree well with 
experiment?0,31,32, 

Some insight into the two types of creeping waves may be garnered from 
consideration of wave propagation along an infinite plane interface between two 
different isotropic media. (Uberall2 contains an indepth review of circumferential 
wave.) If one of the media is a vacuum and the other an elastic solid, Rayleigh? 
showed that an inhomgeneous surface wave results which exponentially decays 
normal to the surface and propagates along the surface. This wave, the classical 
Rayleigh wave, is a combination compressive-shear type wave which propagates 
with a unique velocity, Cp, along the surface. The Rayleigh wave speed varies 
from 0.86 to 0.96 of the shear wave speed in the elastic media as Poisson’s ratio 
varies from 0 to .5. No waves exist in the vacuum. Love34 and Stoneley>> 


extended the analysis to liquid-solid and solid-solid interfaces. The surface waves 
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generated in these situations are often called generalized Rayleigh waves. 
Generalized Rayleigh waves differ from the ideal Rayleigh wave in that the gener- 
alized Rayleigh wave exist in both media and are exponentially damped in the di- 
rection of propagation (nonpropagating waves). The generalized Rayleigh wave 
also decays exponentially in the normal direction from the surface in both media. 
The existence of the wave in the second media results from the seepage of energy 
from the first media which causes the wave to be nonpropagating along the surface. 
(The seepage of energy has resulted in the generalized Rayleigh wave being re- 
ferred to as the Leaky Rayleigh wave.) In addition for the generalized Rayleigh 


waves, new modes of the surface wave exist. The dispersion equation is 
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where Cy is the longitudinal wave speed, Cry is the transverse wave speed and fg is 
the density of the elastic solid, Cr is the propagation speed and pf is the density of 
the liquid (fluid) and C is the generalized Rayleigh wave speed to be determined. 


Equation (1-4) has one real root and the remaining roots are complex. As pr ~> 0, 
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there is one complex root whose real part approaches the Rayleigh wave speed, 
Cr. (The fact that the roots are complex implies that the wave is attenuating in 
the propagation direction.) The real root indicates a propagating surface wave 
which travels mainly in the liquid. The wave is still attenuating normal to the 
surface. The phase speed of the propagating surface wave is less than the 
propagation speed of the liquid (see Viktorov°® for detailed derivation of the roots 
of equation (1-4)). This wave is generally referred to as the Stoneley wave. 
Uberall, et al.39 demonstrates the Watson-Sommerfield transform for the pressure 
release cylinder in a fluid. The transform yields poles which correspond to the 
propagation phase speeds less than the propagation speed of the fluid. These poles 
represent the different modes of the geometric or Franz waves. Uberall, et al. also 
shows that the first Franz wave for the pressure release cylinder is significantly 
more attenuated than the corresponding first Franz wave of the rigid cylinder. 
This attenuation explains why the pressure release backscatter form function does 
not show the classical form function structure associated with the infinite rigid ai 
inder. The analogy is that the superposition of the Franz modes result in the 
unique flat plate Stoneley wave. The attenuation of the Franz wave is a curvature 
effect caused by the tangential radiation of energy into the fluid for each mode. 
Neubauer“? demonstrates the Watson-Sommerfield transform for a solid alumi- 
num cylinder. The transform produces two sets of poles referred to as Franz type 
and R type (Rayleigh type) poles. The R type poles result in dispersive wave 


propagating at phase speeds greater than the propagation speed of the fluid. Frisk 
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and Uberall?’ related the first R type pole, Ry, to the generalized Rayleigh wave 
since its phase speed approaches the Rayleigh wave speed, Cp, as kat». The other 
R type poles are called "whispering gallery” waves and become lateral waves as 
ka>a, 

As stated above, the R-type waves are not present for the pressure release and 
rigid cases. When an incident pulse is used, Neubauer? shows the Franz wave 


and the specular reflection can be time separated with 


Rees e.. (1-5) 
C C 
W p 


where oe is the phase speed of the Franz wave and C,, is the propagation speed 
of the water. As illustrated in figure 1-1, the first term of equation (1-5) 
corresponds to time to travel twice the distance dj in the fluid and the second term 
of equation (1-5) is the time to travel half of the circumference of the cylinder at 
the phase speed of the Franz wave. Due to attenuation, only the first Franz wave 
is of significance (see Neubauer2‘). It is also possible to recover the Franz wave 
by decomposing the scattered pulse into a specular reflection and a creeping wave 
using the assumption that the specular reflection has the same frequency content as 
the incident pulse. This technique is illustrated by Rudgers?8 for a mgid sphere. 
Another possible contributor to the scattered field is multiple internal 
reflections of the shear wave. Neubauer and Dragonette2? demonstrated the 


significance of these multiple reflections. Brill and Uberall*9 also provided 
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theoretical support for the shear wave reflections. 

In order to understand the scattered field of the cylindrical shell, the relation 
between the circumferential waves in the shell and the scattering form function 
must be developed. From the form function versus ka graph for the cylindrical 
Shell, three well defined background regions are apparent. In region 1 or soft 
background region (small ka values), the specular reflection is inverted relative to 
the incident signal and the resonances appear as dips (decreases) in the form 
function on the form function verses ka graph. The second region or transition 
region occurs when the resonances appear as both dips and peaks on the form 
function verses ka graph. In the third or rigid background region, the specular 
reflection is in phase with the incident signal and the resonances appear as peaks 
on the form function verses ka graph. For the shell, Neubauer? states that the 
Franz waves do not measurably affect the form function. The R type waves are not 
present Since they require a semi-infinite medium (e.g. the thickness of the shell is 
usually less than a wavelength and can not support the R type waves). Horton, et 
al.41 made the first attempt to relate motion of the shell to plate theory. Diercks, 
et al.19, established the existence of two groups of circumferential waves in shells 
with different group velocities. They were referred to as longitudinal and flexural 
waves where the group speed of the longitudinal waves is greater than the group 
speed of the flexural waves. Some confusion existed due to the variation in speeds 
that were obtained for the fast (longitudinal) circumferential waves. Herrmann 


and Mirsky42 determined that shell analysis can be divided into two methods. For 


a! 
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very thin shells, the motion is dominated by membrane and curvature effect, 
therefore membrane theory is used. For thicker shells, the motion is dominated by 
flexural, rotatory inertia, and shear-deformation effects, therefore plate theories 
are used. For wavelengths less than A;, the membrane theory is used. 


2 
Ay a (1-6) 


| Ail owe 


where R is the mean radius, h is the shell thickness and vy is poisson’s ratio. For the 
experiments to be conducted in this investigation, the shell is of sufficient thickness 
to require plate theories. The two circumferential groups of waves observed by 
Diercks, et al.19 correspond to symmetric Lamb modes (longitudinal) and 
antisymmetric Lamb modes for infinite plates (wheresymmetry is defined relative 
to the axis located in the center of the plate, parallel to the surface). In order for 
either of the modes to be to be excited, the phases speed of the mode must be 
greater then the propagation speed of the fluid (water). The angle at which the 


mode is excited is determined from 


C 
sin(@.) = a where : is the phase speed of the mode. 


p 


It is possible to generate the first symmetric mode for all frequency thickness (fh = 
shell thickness x frequency) values in water. For a complete analysis, including 


dispersion curves, of mode generation in steel (ratio of longitudinal to shear speeds 





15 


equal to 1.8), see Grigsby and Tajchman‘9. Viktorov2® shows the first symmetric 


Lamb mode phase, Vs.p» and group , Vs,g» speeds vary from 


E 
"s(p.) 7 na oe 


to the Rayleigh wave speed, Cp, as fh varies from 0 to». The first antisymmetric 


mode phase speed, Vap» varies as 


4 
Vue =. = (1-8) 
mel 3p(1 - v) 2 


(approaches 0 as fh-~) and approaches the Rayleigh wave speed, Cp, as fhoa. The 





first antisymmetric mode group speed is equal to twice Vap and also approaches 
the Rayleigh wave speed, Cpr, as fhom. As the frequency thickness increases, 
additional modes of the symmetric and antisymmetric Lamb waves become excit- 
ed. The cutoff frequencies of the symmetric modes are given by the plate thickness 
corresponding to one half wavelength intervals of the longitudinal wave and whole 
wavelength intervals of the transverse wave. The cutoff frequencies of the 
antisymmetric modes are given by the plate thickness corresponding to one half 
wavelength intervals of the transverse wave and whole wavelength intervals of the 
longitudinal wave. The higher modes approach the transverse wave speed of the 


material as the frequency thickness tends to infinity. Viktorov-® states the 
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asymptotic forms of the phase speed for the symmetric and antisymmetric modes of 


order nas 





2 mane 
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where D = (wh)/(2C7). 

The significant observations (discussed in Neubauer?) from the plate theory 
predictions are: 

1. for thin shells, the first symmetric Lamb mode dominates in the low ka 
region and will cause significant dips in the form function verses ka graph at 
intervals of Aka = 3.7. These dips result from the addition of the specular 
Teflection, which is 180 degrees out of phase from the incident signal, and the first 
symmetric mode circumferential wave (at the frequency where the wavelength of 
the circumferential wave corresponds to an integer multiple of the circumference 
of the cylinder), which is in phase with the incident signal, producing a cancelling 
effect in the overall scattered signal. As the shell thickness increases, the dips 
become less pronounced and additional symmetric modes may become excited. 
The excitation of additional modes explains the early confusion in the literature 


involving the phase speed of the fast circumferential waves. 
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2. the first antisymmetric mode is usually not excited in thin shells until larger 
ka values (for aluminum and an inner to outer radius ratio of .98, the first 
antisymmetric mode is not excited until a ka of approximately 55). As the shell 
thickness increases, the frequency required for excitation of the first antisymmetric 
mode decreases (same fh value). (Neubauer shows that for an aluminum shell with 
an inner to outer radius of .9, the antisymmetric mode is observed at a ka=10.) 
The antisymmetric mode circumferential waves cause resonance peaks in the form 
function versus ka graph. The peaks occur because the shell is scattering in a rigid 
manner since the specular reflection is in phase with the incident signal and the 
antisymmetric wave, which is an elastic wave, is also in phase with the incident 
signal. When the wavelength of the antisymmetric mode is the correct value to 
how exact phase matching, the antisymmetric mode circumferential wave and the 
specular reflection add to cause a resonance peak. 

An alternate method of determining the scattered field is to assume a Series 
solution in form of inverse integer powers of ka. This method, called the 
Luneberg-Kline method“4, is adapted from optical theory. The scattered pressure 
iS 


@ 


Vn(r) 
Ps(r) = exp(uky) ) a (1-10) 
(tka)f 


n=0 


Since only integer powers are considered, the circumferential waves are neglected. 
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Due to the difficulty in determining V,(r), only the first two terms are retained, i.e. 
n=2. Keller et al.45 applies this method to the pressure release, rigid and 
impedance boundary conditions. George and Uberall46 compare the Luneberg- 
Kline method with the Kirchhoff method for complex impedance surfaces and 
concluded that the Luneberg-Kline method gives a more accurate measure of the 
scattered field. 

In the previous references with the exception of the approximate methods, 
absorption by the scatterer has been neglected. Schuetz and Neubauer®/ 
demonstrated that there is significant error in the scattered field if absorption 1s 
neglected. Scattering of sound from sound absorbing surfaces has been 
investigated experimentally48-52 and theoretically by using approximate shel] 
theories>>, by assigning a complex impedance boundary condition on the surface of 
the scatterer>4-96 and through the assignment of a viscosity throughout the body of 
the scatter’-8. Since the scatterer to be studied is not homogeneous in the radial 
direction, i.e. layered, the approach by Gaunaurd>? and Flax et al.>8 is most accu- 
rate. This approach solves the exact elasticity equation for the elastic cylinder 
coupled with the exact viscoelastic equation (plane strain case) for the absorptive 
layer. Gaunaurd>® also solves the layer problem using the Sommerteld-Watson 
transform method. The overall result is that the scatterer demonstrates both 
pressure release and rigid-like behaviors at different ka values. After these - 
characteristic behaviors, resonance effects are observed. Gaunaurd’s»’ material 


properties do not correspond to any known materials, so there can be no experi- 





a 


mental verification of his results. Flax and Nuebauer-® solved the same problem 


as Gaunaurd with a lucite coating but, again, there is no experimental verification. 


Numerical Methods 


There are a number of methods which can be used to determine the acoustic 
field scattered from a surface. It is often convenient to reformulate the scattered 
problem as a radiation problem. Chertock®9 proved this can be done with no loss 
of generality. Since most of these methods have been derived for the radiation 
problem, i.e. a vibrating surface in an infinite medium, this literature review will be 
based on this reformulation. 

The radiation problem is conventionally posed as follows (see figure 1-2). 
An arbitrarily shaped finite object with a total surface area S, is contained in an 
infinite, ideal, homogeneous fluid. The density of the fluid, p, and the speed of 
sound, c, in the fluid are known. The region outside the surface is designated Ry 
and the interior region is designated Rj. An arbitrary point in R, is designated x 
and an arbitrary point in Rj is designated x;. Points on the surface S are designated 
X,. The positive unit normal vector ng is directed from a point on S into Rg. Only 
the steady state condition will be considered which implies the pressure and 
velocity are harmonic functions. The harmonic time dependence will be 
Suppressed in this review. The boundary value problem requires that the pressure, 


p(x), satisfy the Helmholtz equation for all x in Ro. 








Figure 1-2 General scattering Geometry for Numerical Methods 
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(V2 + k2)p(x) = 0 (1-11) 


The surface boundary condition is obtained from the normal component of 


Euler’s equation 


ee tg (Vp(x)) a) = -LWP V(x.). (1-12) 
Ss 


The pressure is also required to satisfy the Sommerfield radiation condition at 


infinity 


ee r | EM) - ukp(x) | = 0. (1-13) 


The most straightforward method is to solve the wave equation directly using 
the classical method of separation of variables. This method, however, can be used 
only when the vibrating surface conforms to a separable coordinate system. (See 
Morse and Feshbach®! for examples.) The normal mode method (Rayleigh series) 
for scattering works for those special cases when separable coordinate systems are 
usable for the radiation problem. A disadvantage of the normal mode method is 
that the calculation of the solution functions is not trivial nor always 
straightforward. Approximate solutions to the wave equation have been developed 


but these are extremely situation dependent with little generalization. 
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Due to the difficulty is solving the wave equation directly, alternate integral 
methods have been developed. During the 1960’s four integral equation methods 
came into prominence. The first is the simple source formulation which has its 
Origins in classical potential theory. The basic premise is that the pressure 
(potential) field can be represented by the integral over the surface of a simple 
source density function, o(x,), times the free space Green’s function. It is assumed 


that o(x,) is originally unknown. 


ukr 


p(x) = -Lup | o(x.) - d5(x_) (1-14) 
S 





The integrand of the above integral is singular as x approaches the surface but 
Kellogg52 showed that the limit of the integral does exist and is continuous. The 
normal derivative of the integral is not continuous at the surface and must be 
accounted for. If the surface boundary condition is applied to an arbitrary point x, 


on the surface, equation (9) becomes 





va) + 004) +] | ob) - & | asta, (1-15) 
S 


where the integral is an improper integral. Chen and Schweikert®3 demonstrated a 
method of subdividing the surface into N triangular elements and assuming 0(x,) is 


constant over the each element. This formulation allows the integral to be written 
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in terms of N linear algebraic equations. Chen and Schweikert and Baron et a].64 
work several examples of this method for vibration of the acoustic medium and for 
vibration of elastic solids and shells in contact with the medium. 

The remaining three methods are based on the Helmholtz Integral Method 
(HIM). The Helmholtz Integral determines the pressure at an external field point, 


x, in terms of the surface pressure, p(x,), and the normal surface velocity, v(X<). 


4x p(x) = | | | p(x.) an) g(x,x_) - bwp v(x.) g(X,x,) | dS(x_) 
| S 


where gix,x.) = (1-16) 


The normal surface velocity and surface pressure are related by two integral 
equations called the Surface Helmholtz Integral Equation (SHIE) and the Interior 
Helmholtz Integral Equation (IHIE). A complete discussion of the HIM may be 
found in Baker and Copson®°. Chertock®6 points out a low and high frequency ap- 
proximate relation between the surface pressure and surface velocity. In the low 
frequency limit the in phase component of the surface pressure is considered 


negligible and 
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p(x.) = LW v(x.) m(x.) (1-17) 


where twv(x.) is the local acceleration and m(x,) is the entrained mass per unit 
area. For low frequencies, m(x.) varies slowly and is calculated from incompress- 
ible flow theory. For the high frequency case, the surface is approximated as a 


sphere of radius a, where a is the local mean radius of curvature of the surface. 
p(x.) = ac vix,) [1 - = | (1-18) 


This approximation 1s usually valid for ka>>1. A derivation and additional 
restrictions for this approximation can be found is Morse®/, 

In the 1930’s Kupradze®8 developed what has become known as the Surface 
Helmholtz Integral Equation. To develop this equation, the point x in Ro Is 
allowed to approach the surface S. The results in the following two dimensional 


inhomogeneous Fredholm equation of the second kind 


p(x.) = | | | p(x.) ant) g(x,x_) - Lup v(x.) g(x,x_) | dS(x_) 
S 


where g(x,x_) = (1-19) 


The actual method of solution, i.e. Surface Helmholtz Integral Equation Method 
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(SHIEM), is to solve the surface Helmholtz integral equation for the surface 
pressure given the surface velocity. Then utilizing the surface pressure and surface 
velocity, the HI is solved for the field pressure. Chertock demonstrates a method 
of solution for surfaces of revolution and reduces equation (1-19) to one-dimen- 
sional form. Banaugh and Goldsmith®? describe a related approach for two-di- 
mensional problems. 

The third method pioneered by Copley /0 has become known as the Interior 
Helmholtz Integral Equation Method (IHIEM). It is developed by defining the 


pressure interior to the surface, region R;, as zero for all x; resulting in 


0 = | | p(x.) ant) g(x. 5x.) - LWP v(x.) g(x. 4x.) | dS(x_) 
S 


(1-20) 


This equation, the interior Helmholtz integral equation, defines an acoustic rela- 
tion between the surface pressure and the normal surface velocity which then 
allows usage of the exterior Helmhotz Integral (HI) to determine the exterior field 
pressure. 

The three previous methods have significant difficulties which are described in 
detail by Schenck’1. In summary, the simple source formulation does not always 
exist and the surface Helmholtz Integral equation (SHIE) does not have unique — 
solutions at or near certain characteristics frequencies or wave numbers (these are 


eigenvalues of the interior Dirichlet or Neumann problem). Schenck determines 





26 


additional compatibility relations which must be satisfied for the SHIE to yield a 
unique solution for all frequencies. The Interior Helmholtz Integral Equation 
Method, IHIEM, is unique in principle but exhibits numerical computational 
problems. 

To facilitate usage of the first two methods, knowledge of the characteristic 
frequencies is required so they can be avoided. This requires solving the interior 
eigenvalue problem first. Since numerical approximations are involved in 
discretizing the integrals, a range of values for each characteristic frequency is 
required. The numerical approximation effect becomes increasingly more 
important at higher frequencies due to a reduction in the interval between the 
characteristic wave number associated with the higher frequencies. The overall 
result is that these methods are not generally satisfactory without additional 
modifications or knowledge. 

To eliminate the lack of uniqueness, Schenck’! proposed the Combined 
Helmholtz Integral Equation Formulation (CHIEF). Basically CHIEF discretizes 
the surface Helmholtz integral equation into an N by N system of equations and 
then overdetermines the system by forcing additional compatibility equations to be 
satisfied. The additional compatibility equations are, in fact, the interior 
Helmholtz integral equation applied at selected interior points. Meyer et al./4 
points out two difficulties with CHIEF: (1) there is no simple method to determine 
how many additional interior compatibility equations are required to sufficiently 


overdetermine the system thereby eliminating the nonunique solutions. (2) there is 
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no scheme for choosing the appropriate interior points except to avoid obvious 
symmetry points. Rogers/3 develops a method to determine how many additional 
interior points are required by using eigenvalue decomposition. Recently (August 
1986) Benthein and Barach’4 and Naval Ocean Systems Center, San Diego 
reissued CHIEF in a general program (fortran) format with modern memory 
(storage) allocation procedures which allows rapid program runs with different 
compatibility points so as to arrive at a “best” solution. This new version also 
includes a point source check feature. 

A variation of CHIEF has been proposed by Shaw’. In this method 
additional points are placed in the exterior region, Ro, to overdetermine the 
exterior Helmholtz integral equation. Unlike the interior problem, the potential 
(pressure) at the additional field points is unknown and must be solved by itera- 
tion. 

One method of eliminating the uniqueness problem is to use a different 
Green’s Function. Ursell’® provides a new function which involves the 
computation of infinite series. The shortcomings of this method are primarily com- 
puter related. 

Burton? / points out that a modified differentiated form of the surface 
Helmholtz Integral Equation (SHIE) also describes the surface pressure. It also 
has nonunique solutions corresponding to internal eigenvalues of the associated 
interior equation, but these are different from the eigenvalues of the previous 


interior Helmholtz Integral Equation (IHIE). The differentiated form results from 





28 


taking the normal gradient of the surface Helmholtz Integral Equation (SHIE). 


This yields a Fredholm equation of the first kind. 





2 
Ss 
S 


NS an(x) an(x_) a dae n(x) 
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The significant observation by Burton and Miller /8 was that combining the surface 
Helmholtz Integral Equation (SHIE) and its differentiated form results in a 
solution that is common to both equations. Since the eigenvalues are different for 
each equation, the solution of the combined equation will be unique. This 
combined equation is the starting point for Meyer et al. 72 application to radiation. 
The combined equation is 


d p(x) 
d n(x) 





p(x) + a 


] a 
oe | | p(x.) an(x_) g(x. 5X.) - LWP AC) g(x, xX.) dS(x_) 
S 


] a°9(x,x.) 0 g{x,Xx.) 
ae Qn Des) an(x) dn{x.) : es) d n(x) hae 
S 


(eee) 


where @ is a coupling constant. Meyer et al. points out that the term 
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a°9(x,x.) 
p(x.) an(x) an(x_) dS(x_) | le) 
S 


is strongly singular and can not be numerically integrated. A method of overcoming 
this problem is then presented. Terai’? has also treated the highly singular 
integral by reducing it to a contour integral which is nonsingular for most 
practically realizable surface elements. The selection of a is recommended to be 
-t/k by both Reut80 and Meyer et al. 

In an effort to evaluate which method to use, Tobocman®! conducted a Study 
of scattering from a rigid sphere in the intermediate wavelength region (i.e. 
wavelength comparable to target dimensions) using both CHIEF and the Burton- 
Miller method. The study was performed at four characteristic wave numbers 
between 3 and 10. The conclusion is that CHIEF converges more rapidly than the 
Burton-Miller method but at higher wave numbers, k>10, CHIEF is more time 
consuming and becomes impractically slow. In the frequency range of interest, 
CHIEF out performed the Burton-Miller method both in time considerations and 
accuracy. A Surprising result is that at higher characteristic wave numbers, the 
surface Helmholtz Integral Method (SHIEM) (without overdetermination) showed 
an increase in accuracy and was on the order of the Burton-Miller method. 
Therefore, if precise accuracy is not required at higher wave numbers, the SHIEM 
will give rapid results. Since this was shown for the rigid sphere only, extension to 


other geometries is cautioned. 
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The Simplified-Helmholtz-Integral Program (SHIP) developed by Rogers®2 
solves the Helmholtz Integral (using the SHIEM) specifically for finite circular 
cylinders and free-flooded ring transducers. SHIP takes full advantage of the axial 
symmetry which results in the high speed of this method relative to others. The 
standard approach of the previously discussed methods is to solve the SHIE by 
dividing the surface into small surface subdivisions and solving the surface 
integrals. The axial symmetry allows: 1. the surface to be divided into bands (no 
division 1n @ direction) which reduces the number of elements. 2. the conversion of 
the surface integrals into ordinary line integrals. For application to the scattering 
problem, SHIP requires complete axial symmetry and is, therefore, only usable for 
axial or end on incidence. SHIP is easily adaptable to the pressure release bounda- 
ry condition. 

The follow-on step to the overall scattering problem is to match the acoustic 
radiation condition with the structure reaction or fluid-structure coupling. This is 
accomplished by matching the normal surface velocity defined by the elastic baer 
librium relations with that defined by the equations of the fluid field (HI). 
Mathews®3 developed a program called ADRAD which uses the finite element 
technique for the structural material and the Burton-Miller method for the fluid. 
He strongly recommends that CHIEF be avoided due to the previously mentioned 


problems. 
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Recently little has been done to improve the theoretical treatment of the 
radiation problem although gains have been made in the numerical implementa- 
tion of the integral equations. Schuster and Smith®4and Seybert et al.85 have 
developed a more sophisticated boundary integral equation method of solution by 
using elements that are curvilinear and by representing acoustic variables with 
higher order interpolation functions in lieu of the planar elements and piecewise 
constant acoustic variables. The result is greater accuracy for the same number of 
elements, less computing time, but a more complicated formulation is required. 

Waterman®® developed an alternative method to the HIEM called the T- 
Matrix method (TMM). In this method the potential field is expanded in terms of 
the spherical wave functions which is then substituted into the Helmholtz equation. 
The result is a set of linear algebraic equations for the expansion coefficients. 
Tobocman8/ compares the two methods for a rigid prolate spheroid of various 
aspect ratios. His conclusions are summarized as follows: (1) The TMM 
converges more rapidly than the HIEM for spherical targets. (2) The rate of con- 
vergence deteriorates very rapidly for the TMM with increasing aspect ratios but 
the convergence of the HIEM shows little change. (3) The TMM was severely af- 
fected by overflow and roundoff error. (4) As the number of partial waves is 
increased in the TMM, the matrix which must be inverted becomes increasingly ill 
conditioned which is exactly opposite from the HIEM. Numrich et al.88 has ex- - 
tended the TMM to finite cylinders using cylindrical basis functions. Additional 


results and usages of the TMM may be found in Varadan and Varadan®, 
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An alternative to the previously discussed methods is to utilize variational 
formulations to solve the Helmholtz integral equation. Until the recent work by 
Wu, et al.99 and Pierce?! very little interest in this approach has occurred since the 
1950’s. Morse and Feshbach®! contains a summary of this previous work. Wu, et 
al., demonstrates a new approach of the variational principle to the Helmholtz 
integral equation (i.e. formulation contains all self-adjoint kernels) to transversely 
oscillating disks. Pierce derives a stationary expression for radiated acoustic power 
and target strength. In addition, Pierce discusses some alternatives or tricks to 
enhance the solution of the variational formulation but suggest additional study is 
required for comparison with other computational algorithms. 

CHIEF has proven to be reliable for the finite cylinder particularly since the 
characteristic frequencies are well known. The particular advantages of the 
Benthein and Barach CHIEF program/4 is that it is pre-programmed for numer- 
Ous coordinate systems which allows the program to be applied to an almost unlim- 
ited number of shapes. As the shape becomes more complex, the determination of 
the number and locations of the eigenvalues becomes almost impossible. The 
aforementioned program allows the easy input of numerous interior points which 
partially mitigates this problem. For scattering from the finite cylinder with axial 


incidence, the most accurate and computationally efficient alternative is SHIP®2. 
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Determination of Dynamic Viscoelastic Moduli 


The closed cell rubber materials are considered to be viscoelastic with high 
losses. AS with most rubbers, the dynamic behavior % strongly frequency, loading, 
and temperature dependent. The manufacturing process is also important, in 
particular, the addition of fillers such as carbon black?2 and the cure times. There 
are four distinct classes of methods employed to determine the dynamic viscoelas- 
tic constants: 

1. free vibration method. 

2. direct observation of stress-strain curves. 

3. resonance methods. 

4. wave propagation methods. 
The free vibration method has been extensively utilized since its main advantage is 
simplicity. The method is particularly suitable for measurements at low 
frequencies and for low loss materials. Lethersich?? has demonstrated the free vi- 
bration method for frequencies up to 1000 Hertz using photographic recording 
methods. The direct stress-strain method requires special equipment and involves 
considerable difficulties in determining the stress-strain curve at high loading rates. 
Kolsky?4 points out the advantage of this method is that no assumptions are made 
regarding the nature of the dissipative forces and the linearity of the system. The 


three other methods assume linearity (mechanical behavior is independent of 
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amplitude) and that the dissipative force is proportional to the time rate of change 
of the strain. 

In the audio frequency range, the dynamic moduli and loss factors are best 
determined by resonance methods which are described by Read and Dean? and 
by wave propagation methods. There are basically three types of vibrations 
employed in resonance studies which in order of increasing frequency range, are 
(1) flexural, (2) torsional and (3) longitudinal. For isotropic materials, Read and 
Dean state that the complex Young’s modulus is customarily determined from flex- 
ural —— tests in the frequency range of 20 Hz to 10 kHz and by longitudinal 
resonance tests in the frequency range of 5 kHz to 50 kHz. The torsional reso- 
nance tests give the complex shear modulus for frequencies greater than 500 Hz. 
All of these tests are inherently low amplitude with average strain amplitudes on 
the order of 10°. This low amplitude generally eliminates any problems with 
nonlinearity of the dynamic properties. 

The method of exciting the flexural waves (i.e. free-free, clamped-free etc.) is 
important since significant errors and spurious resonances may be introduced. The 
flexural resonance method originated with Ide?® and has undergone numerous 
improvements as test equipment and theory have evolved. Spinner and Tefft?’ 
provide and excellent summary of the methods and theory. Tefft and Spinner?8 
also provide a review of torsional measurements for rectangular and cylindrical 


bars. 
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The longitudinal resonance method was first utilized by Quimby?? who 
compared the resonances in a viscoelastic bar, determined by Stokes’ theory, with 
the experimental results for glass, copper and aluminum bars. Quimby excited the 
specimen bar by attaching an oscillating quartz (piezoelectric) crystal to one end 
and then measured the particle velocity at the other end. The correlation between 
the experimental and theoretical results was considered excellent. 

Ballou and Silverman!00 modified the technique by using a light source and 
photoelectric cell at the free end to record the displacement. Norris and 
Young!01,102 used the longitudinal resonance method and attached 
accelerometers to both ends. The fixed end was driven by a constant frequency 
source. Madigosky and Lee!03improved on Norris and Young’s technique by us- 
ing miniature accelerometers (Bruel & Kjaer type 8307) and a random noise 
source as the forcing mechanism. The accelerometer signals were analyzed using a 
Fast Fourier Transform (FFT) spectrum analyzer. Both Norris and Young and 
Madigosky and Lee utilized Lee’s!04 observation that if the frequency is adjusted 
such that the two ends are out of phase by 90 degrees, the resonances are easily 
determined. Madigosky and Lee also demonstrate wide applicability to both soft 
(lossy) polymers and rigid metals. Capps!05 has used this method to determine the 
dynamic Young’s modulus of some commercially available polyurethanes. No 
published data for closed cell rubbers or corprene appears to be available 106,107 

Cramer108 compares the Young’ modulus determined by the flexural 


resonance method and by the longitudinal resonance method for rubber rods in the 
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1kHz to 8kHz region. He concluded that they agreed within the experimental 
accuracy. 

Nolle!09 has investigated the elastic properties of rubber-like materials for 
frequencies between 0.1 Hertz and 120,000 Hertz using’a combination of five ex- 
perimental methods (one free oscillation method, three resonance methods and 
one wave propagation method). 

The wave propagation methods have been extensively used for non-dispersive 
materials and ultrasonic measurements. Kolsky94 points out the advantages of this 
method over the other methods are: 

1. a large range of frequencies can be covered with a single specimen 
(specimen sizing is not particularly critical). 

2. in determining the attenuation, losses caused by extraneous supports 
can easily be removed. 

3. the method is highly accurate for non-dispersive materials. 
Two disadvantages of the wave propagation method are: 

1. generation of the particular wave of interest can be complicated. 

2. interpretation of the results for dispersive materials can be very in- 
volved. 
This investigation wil! use wave propagation methods. 

Master curves (complex modulus-frequency and loss factor-frequency curves 
at constant temperature) can be generated for modulus measurements at different 


frequencies and temperature. This procedure is based on the time-temperature 
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superposition principle outlined by Ferry!10. This principle 1s based on the 
viscoelastic material becoming more elastic due to greater molecular mobility as its 
temperature is increased at constant frequency. The converse is also observed. 
The result is that a decrease in temperature at constant frequency causes the same 
change in dynamic properties as an increase in frequency at constant temperature. 
Although there is no limit to the range this data may be shifted, Madigosky and 
Lee suggest that reliable interpolation is assured two decades above and below the 
measured range. 

For rubber-like materials it is well established that Poisson’s ratio, v, 1s 
approximately 0.49. Since the materials to be studied have large concentrations of 
air, it is expected that the Poisson’s ratio will be less than 0.49 but it 1s unclear how 
much less. 

The plane wave modulus, M, can be determined from the longitudinal wave 
speed. The common method of determining the wave speed is by using pulse and 
comparing arrival times with and without the material present. The attenuation, 
corrected for spreading loss and reflection, will yield the complex wave speed in 
the material. 

The reflection and transmission coefficients are experimentally determined 
from panel measurements at various angles of incident (panel measurements 
involve placing a flat panel of the material with a known thickness in a sound field - 


and measuring the transmission through or reflection off the panel. The panel 
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dimensions are usually large enough so that the panel may be considered one 
dimensional, i.e. an infinite planar surface). The experimental technique is 
documented by Barnard et al.111 and the theary is based on the work of 
Brekhovskikh!12, Folds and Loggins!13 corrected the errors in Brekhovskikh and 
Barnard et al. and extended the theory to multiple layers using a matrix approach. 


To predict the elastic behavior of the material, two moduli must be 


_ determined. The plane wave modulus is determined from panel measurements 


and the Young’s modulus, E, is determined from wave propagation in a bar. 
Losses in the material are accounted for by allowing the elastic constants to be 
complex (now, they are technically viscoelastic constants). The following 
assumptions concerning the materials are made: 

1. The material 1s isotropic. 

2. The air inclusions are relatively uniform in size and dispersion 
(homogeneous). 

3. The air inclusions are small relative to a wavelength. 

4. The material is stressed within the elastic limit. 

5. The amplitude of the acoustic displacement is small enough so as to 
not affect the dynamic moduli. 
The first assumption implies that the elastic properties are the same in all 
directions. Homogeneous implies that the results will be sample independent. 


The third assumption allows the material to be treated as a composite. The fourth 





— 


By 


assumption allows the material to be modeled with classical elasticity theory. The 
final assumption eliminates displacement amplitude as an independent variable n 


the determination of the elastic constants for small amplitudes! 14. 
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CHAPTER II 


Experimental Procedures 


The experimental procedures are divided into two categories: determination 


_ of viscoelastic constant and acoustic Scattering experiments. 


——_——— = —<.-----_= = 


Determination of Material Viscoelastic Constants 


The materials investigated are grouped into two classes: closed-cell rubbers 


and rubber-cork composites. All materials are assumed to be isotropic. The four 


materials are: 


neoprene rubber with air cells 


(RUBATEX™ No. R-411-N) 


nitrile rubber with air cells Os cCace i muapers 


(RUBATEX™ NO. R-416-N) 


neoprene rubber with cork 
(GROENDYKE~ No. 712 [corprene]) 


nitrile rubber with cork ODA SCOTS (Seln[ AS oe 


(GROENDYKE*™ No. 725) 


(" refers to RUBATEX Corporation and GROENDYKE Manufacturing 


Company, corprene is the tradename of the ARMSTRONG Rubber Company) 
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In order to fully characterize the material, two viscoelastic constants and the 
density are required. The viscoelastic constants will be treated as complex elastic 
constants. The density is determined by the standard method of weighing a known 
volume of the material. The measured densities are all within the manufacture’s 
specified ranges. The constants that are determined are the complex Young’s 


modulus, E, and the complex plane wave modulus, M. 


Complex Young’s Modulus 


The complex Young’s modulus is determined from the relation 


2 
E= pc (2-1) 


where Cz is the complex bar speed. The experimental setup is shown in figure 2-1. 
A long rectangular strip of the material is glued to a Bruel and Kjaer (B&K) 
shaker which is driven by a broadband pulse produced by a Wavetek model 275 
signal generator and amplified by a Kronhite Model 7500 amplifier. The lateral 
dimensions of the strip are approximately .25 inches with a length of greater than 
30 inches. Since the length is very much larger than the other two dimensions, the 
Significant propagating wave is the extensional wave which is governed by the 


Young’s modulus. 
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The procedure is to record the broadband pulse waveform at three locations 
along the sample and at several center frequencies. The filter settings remained 
constant for each center frequency group so as to ensure a consistent time delay. 
Each waveform was time shifted to place the earliest pulse for a particular center , 
frequency group near the time origin. This shift is somewhat arbitrary and not 
always the same for all groups. It is of paramount importance that the time shift be 
the same for all samples within a particular center frequency group. A fast Fourier 
transform (FFT) (standard IMSL FFT routine used) is then performed on the 
shifted waveform. The purpose of the time shift is to reduce the frequency of 
discontinuities in the FFT phase caused by the signal being time displaced from the 
origin. Although it is not theoretically required, the reduction of the numerous 
phase discontinuities facilitates unwrapping (removal of the 2a discontinuities 


resulting from the tangent function). The phase speed is determined from 


wl | 

Cig + nae 2-1) 
where Ag is the unwrapped phase difference and L is the distance between two 
measurement points. Since the exact multiple of cycles between the two sample 
points is not obvious, the n2x multiple is required for proper scaling. The 22 
multiple is determined from the time domain waveform. The time difference 


between two apparently similar points on the two pulses is measured to obtain an 


estimate of the phase speed. Ideally the center of the pulse should be used for 





estimation, however this point is often difficult to ascertain; therefore the largest 
peak was used and assumed to be at the center frequency. An error in the 
selection of n is easily recognized since there is a significant deviation in the calcu- 
lated phase speed from the estimated value. The distance between sample points 
was kept small so as to prevent significant distortion of the waveform caused by 


dispersion. Specifically, the 2% multiple is determined from 





n2n = = - Ag (2-2) 





Amp 
] 2 
@=- a Yog | | (2-3) 


The FFT amplitudes also determine the boundaries of the usable frequency range 
for the individual center frequencies. There is no specific guidance as to what 
amplitude decrease from the maximum results in meaningless data. The 
procedure employed is to use data with amplitudes greater than 10% of the 
maximum (center frequency) amplitude and evaluate the data using a 
"reasonableness criteria". The "reasonableness criteria" is based on: 

j, there is sufficient overlap between frequency bands that the values should 
be approximately the same. 


2. the speed and attenuations of these materials are slowly varying 
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(increasing) with frequency. 

3. the frequency spacing is small enough so as to ensure large 
variations between individual frequency points do not occur. 

The wave was recorded at each sample location using a stereo phonograph 
needle as a displacement sensor. The two outputs from the needle were summed 
using an Ithaco preamplifier operating in the differential input mode which 
maximizes the normal displacement. The amplitude of the signal is not dependent 
on the contact pressure provided contact is maintained and the pressure is not 
great enough to cause the needle to contact the backing assembly which results in 
distortion of the observed waveform. 

The sample was suspended vertically with the shaker at the top. The bottom 
of the sample was glued to a 2 inch square lead block simulating a rigid boundary 
condition. Attenuation in all the test materials was large enough to require the 
sample points be located near the shaker. The attenuation and large distance from 
the rigid end allows multiple reflections to be neglected. 

Since temperature control facilities were not available, complete master 
curves were not obtained. Extrapolation of the data using a log-log plot has been 
shown to be valid for up to 2 decades above and below the measured data (see 
Madigoski and Lee103), The data was obtained over the approximate frequency 
range of 300 Hz to 3000 Hz. The upper frequency is limited by attenuation causing 
the observed pulse to be dominated by noise. The maximum frequency required 


for the shell modeling is 60,000 Hz which can be obtained by extrapolation. 
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The above measurement procedure is subject to the following sources of 
errors: 

1. errors in digitizing the data. 

2. computer round off error during the FFT procedure. 

3. misalignment of the sample over the center of the shaker causing introduc- 
tion of surface and flexural waves. 

4. error in determining the distance between sample points resulting from 
uncertain positioning of the stereo needle. 
Errors introduced by means 1 and 2 are random and assumed to be very small 
resulting in a less than 1% combined error. The generation of surface and flexural 
waves can be disregarded since they travel slower than the extensional wave and 
are highly attenuated. It is estimated that the sample point locations are accurate 
to with in 2mm resulting in a possible 10% error for small sample point separations 


at the higher frequencies. 


Plane Wave Modulus Determination 


The plane wave modulus is also determined by traveling wave techniques. 

The experimental apparatus is illustrated in figure 2-2. A broadband pulse is 
generated by a USRD F-33 transducer and measured by a B&K 8103 hydrophone. 
The 8103 hydrophone is placed on the central axis of the F-33 transducer. A 36 - 
inch by 36 inch panel of .25 inch thick sample material is placed between the F-33 


transducer and the 8103 hydrophone. The sample is attached to a wood frame for 
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support and is weighted so as to suspend vertically. The panel is positioned 15 
inches (38 cm) from the F-33 transducer. This places the panel in the far field of 
the transducer for all measured frequencies. Far field measurements are utilized 
SO as to avoid the axial (z axis) field amplitude variations that occur in the near 
field of a piston-type transducer such as the F-33. In addition, far field 
measurements allow the usage of a receiving hydrophone (8103) which is physically 
smaller than the projector (F-33 transducer) since the x-y field variations that exist 
in the near field are mitigated. For a frequency of 50 kHz, transition to the far 
field occurs at 13.1 inches (= [F-33 radius]*/)). As illustrated in figure 2-3, for 
15kHz the columnar portion of the axial field (near field) extends for a distance of 
3.94 inches (10 cm) from the F-33 transducer. This results in a spreading angle of 
45° from the center axis; therefore at the panel location of 15 inches from the 
transducer, the sound will cover 15 inches of the panel in the x-y plane. Diffraction 
by the edges of the panel could interfere with the received signal for frequencies 
less than 15kHz. | 

The experimental procedure is to record the waveform from the 8103 
hydrophone with and without the panel in place. An FFT is performed on the 
signal as was done in the Young’s modulus experiment. From a comparison of the 
FFT data, the phase speed and attenuation is obtained. The closed cell rubbers 
are highly attenuating and multiple internal reflections within the panel can be 


neglected. The phase speed, c,, is calculated from 





49 






F-33 Transducer 





10 cm 
Panel 
Near Field 


Figure 2-3 F-33 Beam Spreading at 15 kHz for Panel Measurement 
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_— (2-4) 


where L is the panel thickness (0.25 inches (0.635 cm)), cy is the phase speed in 
water (1.5 x 10° cm/sec), and Ag isthe unwrapped phase difference between the 


two signals. The attenuation, a, is determined from 


gah panel : mike -aL 
<a oS el 
without panel eZ 
where ry is the specific acoustic impedance of the water (= p,c,,)and ro is the 
specific acoustic impedance of the panel (p7c7). The internal multiple reflections 
for the corprene and cork-nitrile composites can not be neglected and are removed 
by an iterative method described in Appendix E. 
There are four possible sources of error in this experimental procedure: 
1. errors in the digitizing of the data. 
2. computer round off error during the FFT procedure. 
3. misalignment of the F-33 transducer, panel,or 8103 hydrophone. 
4. diffraction effects. 
The first two error are the same as for the Young’s modulus experiment and are 
considered very small. Misalignment resulting from the F-33 transducer and the 


panel not being exactly parallel results in a variation of the effective thickness. A 


10° misalignment will cause an error of 1.5% in the speed measurement. 
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Additionally, shear and flexural waves can now be generated in the panel. The 
fourth source of error, diffraction effects, is avoided by setting the lowest usable 
frequency greater than 15kHz. Since the wave generation effects will not be 
investigated, the maximum error assumed will be equal to the maximum data 
spread for the same frequency between two adjacent center frequency data groups. 


The maximum observed data spread in the phase speed measurement was 7%. 


Scattering Experiments 


Two basic scattering experiments were performed: (1) scattering at normal 
incidence from a long circular tube simulating an infinite cylindrica] shell, and (2) 
Scattering from a thick, finite cylindrical shell with flat end caps. The pertinent 


distances for the scattering experiments are determined acoustically. 


Infinite Cylinder Experiment 


For the infinite cylinder simulation, the tube is commercial grade, cease 
304 stainless steel tubing with an outside radius of 1.5 inches (3.81 cm) and a wall 
thickness of 0.035 inches (0.089 cm). The ratio of inner to outer radii is 0.977 
which implies the tube should scatter like a thin shell. The length of the tube is 8 
feet. 

The experimental setup is illustrated in figure 2-4. The acoustic source is a 
USRD F-33 transducer. The receiving hydrophone is a B&K 8103 hydrophone. 


Both the cylinder and F-33 are suspended vertically which ensures they are 





F-33 Transducer 8103 Hyrophone 


56 cm 90.3 cm 


dj = 56cm + 90.3 cm = 146.3 cm 


d> = 90.3 cm 


Figure 2-4 Experimental Setup for Infinite Shell Measurements 
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parallel. The B&K 8103 and target are both located in the far field of the 
transducer. Due to the length of the target and location of the transducer, end 
effects of the target are negligible for center frequencies greater than 18 KHz. 
Since broadband pulses were used, the scattered signals from the ends could be 
time-gated out for center frequencies less the 18 KHz. The measurements are tak- 
en over an approximate ka range of 1.5 to 15 using ten broadband pulses. The 
reference (incident) signal is obtained by replacing the target with the B&K 8103 
hydrophone (at the same physical location). ~The FFT amplitude of the reference 
signal is compared to the FFT amplitude of the scattered signal to determine the 
form function. 

The scattered signal must be corrected for spherical spreading in the vertical 
plane since the incident wave is spherical not planar. The form function corrects 
for cylindrical spreading in the horizontal plane. Spherical spreading in the 
backscatter direction can be represented as cylindrical spreading in two orthogonal 


planes; therefore the scattered signal is multiplied by 


correction = ———S == (2-6) 


where dj] is the distance from the F-33 transducer to the center of the target and d> 
is the distance from the center of the target to the 8103 hydrophone. 


The following four experiments are performed: 
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1. uncoated shell only. 

2. .25 inch (.635 cm) nitrile coated shell. 

3. .0625 inch (.159 cm) corprene coated shell. . 

4. .3125 inch (.794) corprene coated shell. 

The coatings were glued with contact cement and care was taken to prevent any air 
inclusions between the coating and the shell. The seams of the coating were placed 
on the forward scatter side of the target. 

The measurements on the uncoated shell were used to validate the experi- 
mental method and theoretical calculations since similar results are well 
documented. The purpose of the two thicknesses of corprene was to determine if a 
pressure release boundary condition could be achieved and the effect of thickness 
towards the realization of this goal. A secondary purpose of the two thicknesses of 
corprene was to validate the theory for low attenuation material and for thickness 
variation. The nitrile coating’s main purpose was to attempt to achieve a pressure 
release boundary condition and secondary purpose was to validate the theory for 


large attenuation materials. 


Finite Cylinder Experiment 


The finite cylinder is constructed of commercial grade, seamless, 304 
Stainless steel pipe with an outer radius of 2.0 inches (5.08 cm) and a wall thickness 
of 0.1875 inches (0.476 cm). The ratio of inner to outer radii is 0.906. The length 


is 12.75 inches (32.39 cm). The end caps are constructed of 304 stainless steel plate 
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and are 0.75 inches (1.905 cm) thick. The scattering geometries used are: 

(1) incident wave traveling along the center axis (end on), illustrated in 
figure 2-5. 

(2) incident wave traveling normal to the center axis (side incidence), 
illustrated in figure 2-6. 

This cylinder should scatter in a similar manner to the classical thick shell, i.e. the 
first antisymmetric and symmetric Lamb modes in addition to other resonance 
Structure should be present. 

The purpose of this experiment is to determine if the scattered pressure of 
the finite thick shell would as a pressure release finite cylinder as determined by 
the SHIP and CHIEF radiation programs. This requires the characteristic 
scattering behavior of the stainless steel structure to be eliminated by coating the 
finite cylinder with a highly attenuating, low acoustic impedance coating. In other 
words, can the structure be made to "look" pressure released thereby acoustically 
decoupling the structure by use of the coating. 

For each geometry, measurements are taken for the uncoated cylinder and 
for the cylinder coated with 0.25 inches (0.635 cm) of the closed cell neoprene. 
The ka range covered is approximately 1.5 to 15 by using broadband pulses as in 
the previous experiments. An FFT is taken on the scattered signal and the FFT 
amplitude is compared to a common reference obtained at the geometric center of 
the cylinder. Positioning information is not of sufficient accuracy to make the FFT 


phase information useful. For both geometries, the 8103 hydrophone is located in 
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F-33 Transducer 8103 Hyrophone Shel] 


56 cm 90.3 cm 


R = Radius of finte cylinder (5.08 cm) 


Figure 2-5 Finite Cylinder Experimental Setup for Axial Incidence 
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F-33 Transducer 8103 Hyrophone | Shel] 


56 cm 90.3 cm 


R = radius of finite cylinder 


Figure 2-6 Finite Cylinder Experimental Setup for Normal Incidence 
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the far field of the target. 

Experimental errors arise predominantly from misalignment of the target. 
Two possible errors are: 

(1) the target is not positioned at exactly the same location relative the 
reference location. 

(2) the target is angularly misaligned so that the incident signal is not exactly 


norma! to the surface of interest. 





>) 


CHAPTER ITI 


Experimental Results and Discussion 


Viscoelastic Moduli Determination 


For the Young’s modulus and plane wave modulus determinations, the 
measured quantities are the phase speed, Cp and attenuation, a. These two 
quantities combine to give the wave number of the particular wave, either exten- 


sional or longitudinal wave, generated. 


k = = + Ula (3-1) 
p 


The attenuation is determined in units of nepers/cm and converted to dB/cm (by 
multiplying by the conversion factor of 8.6859). The respective viscoelastic 


constant (complex quantity) is determined from the complex wave number by 


ao (3-2) 





The measurements were not performed on numerous samples or lot numbers so 


direct application of these results to other samples is cautioned. Accuracy of the 
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results is difficult to estimate. However, the data spread obtained for different 
center frequency groups provides a reasonable measure of the accuracy. The 
Kramers-Kronig relation, linking the attenuation with dispersion (frequency 
dependence of the phase speed), is utilized as a check of the data for a particular 
measurement sequence. O’Donnell, et al.115, derives the local Kramers-Kronig 


relation for linear acoustic systems: 


= (oes) 


where f is the frequency, C, is the phase speed, and dC, /df is the experimentally 
measured local slope from the phase speed verses frequency curve at Cp. The 


densities of the materials are tabulated in table 3-1. 


Table 3-1 


Density of Materials 


Material Density (grams/cc) 
closed cell neoprene 0.224 
closed cell nitrile 0.248 
corprene 0.720 


cork-nitrile 0.670 
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The closed cell neoprene data is shown in figures 3-1 through 3-4. 

The extensional data becomes unusable at frequencies greater then 3500 bizare 
maximum spread of data for the extensional speed measurements was 10% and for 
the extensional attenuation measurements was less than .5 dB/cm. The maximum 
spread of data for the plane wave speed measurements was 5% and for the plane 
wave attenuation measurements was less than 1.0 dB/cm. The extensional 
attenuation predicted by the Kramers-Kronig relation using the experimentally 
determined phase speed - frequency curve is within 7.5% of the experimentally 
measured attenuation. The plane wave attenuation predicted by the Kramers- 
Kronig relation is within 3% of the measured attenuation. Extrapolation of the 
extensional data to higher frequencies indicates that the Poisson’s ratio is slightly 
greater than 0. Microscopic inspection of a cross section of the material indicated 
the material should be isotropic. 

The closed cell nitrile data is shown in figures 3-5 through 3-8. The 
attenuation is too large to allow accurate measurement of extensional data at 
frequencies greater then 3500 Hz. The maximum spread of data for the extension- 
al speed measurements was 5% and for the extensional attenuation measurements, 
the maximum spread of data was less than .5 dB/cm. The maximum spread of data 
for the plane wave speed measurements was 4% and for the plane wave attenua- 
tion measurements was 1.0 dB/cm. The extensional attenuation predicted by the 
Kramers-Kronig relation is within 3.7% of the experimentally measured 


attenuation. The plane wave attenuation predicted by the Kramers-Kronig rela- 
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tion is within 1% of the measured attenuation. Extrapolation of the extensional 
data to higher frequencies indicates that the Poisson’s ratio is between 0 and 0.2. 
Microscopic inspection of a cross section of the a .25 inch thick material indicates 
the material is isotropic. Microscopic inspection of the 1/16 inch thick nitrile 
reveals a very different air bubble structure than is present in the .25 inch nitrile. 
The thicker nitrile consists of relatively small, uniformly dispersed air bubbles 
throughout the matrix whereas the 1/16 inch nitrile has the same small bubble 
consistency but with the additional presence of the very large bubbles randomly 
dispersed throughout the matrix. It is very questionable if the measurements for 
the .25 inch thickness apply to the thinner material. 

The corprene data is shown in figures 3-9 through 3-12. The extensional 
attenuation is too large at frequencies greater than 5000 Hz. to allow accurate 
measurement. The maximum spread of data for the extensional speed 
measurements was 3% and for the extensional attenuation measurements was 0.7 
dB/cm. The maximum spread of data for the plane wave speed measurements was 
1.5% and for the plane wave attenuation measurements was less than 0.3 dB/cm. 
The extensional attenuation predicted by the Kramers-Kronig relation is within 
25% of the experimentally measured attenuation. The Kramers-Kronig relation is 
not easily applied to the plane wave data since the slope of the phase speed-fre- 
quency curve is small and not obvious. Small variations in the slope result in 
differences of 50% of the predicted attenuation. As a check of the data, the FFT 


analysis was performed on the early portion of the transmitted waveforms where 





40000 


SPEED (crr/sec) 
30000 32500 35000 37500 


27500 


0 


CORPRENE EXTENSIONAL SPEED 


1000 2000 3000 4.000 5000 
FREQUENCY (Hz) 


Figure 3-9 Corprene Extensional Speed versus Frequency 
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. Figure 3-10 Corprene Extensional Attenuation versus Frequency 
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the length of the waveform was selected so as to prevent any multiple internal 
reflections from interfering. The above analysis of the shortened waveform was in 
excellent agreement with the analysis of the complete waveform which accounted 
for multiple internal reflections. Horsley and Thompson!16 report the longitudi- 
nal speed of corprene (manufactured by Armstrong Rubber Company) to be 
3.0x10° cm/sec at 10 kHz as determined by impedance tube measurements. A\I- 
though these two materials are not produced by the same manufacturer, they do 
have the same military specification (mil spec). The extensional phase speed is sig- 
nificantly greater than the longitudinal phases speed for comparable frequencies 
(well beyond any possible experimental error) which indicated the Young’s modu- 
lus is larger than the plane wave modulus. The moduli are related by Poisson’s ra- 
tio as follows 


Eh =) 


(ieee te 2)) ae 


For any positive value of Poisson’s ratio, the plane wave modulus is always greater 
than the Young’s modulus and they are equal at a Poisson’s ratio of 0. The 
conclusion is that the material is slightly anisotropic. Microscopic inspection of 
the longitudinal wave path and the extensional wave path shows a preferential 
alignment of the cork particles so that the thinner portion of the cork particle is 
aligned in the extensional wave propagation direction. This particle alignment 


mostlikely occur-red during the rolling phase of fabrication. Since the plane wave 
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data is directly applicable to the frequency range of interest, a Poisson’s ratio of 0 
is assumed and the Young’s modulus is set equal to the plane wave modulus. The 
shear modulus, which is required for the viscoelastic layer program, is determined 


from 


——— (3-5) 


The validity of this assumption will be determined when the actual scattering 
measurements are compared with theoretical calculations for a Poisson’s ratio of 0 
and 0.4. 

The cork-nitrile data is shown in figures 3-13 through 3-16. The extensional 
attenuation is too large at frequencies greater than 5000 Hz. to allow accurate 
measurement. The maximum spread of data for the extensional speed 
measurements was 3% and the maximum spread of data for the extensional attenu- 
ation measurements was 0.2 dB/cm. The maximum spread of data for the plane 
wave speed measurements was 8.3% and for the plane wave attenuation 
measurements was less than 1.5 dB/cm. The extensional attenuation predicted by 
the Kramers-Kronig relation does not agree with the experimentally measured at- 
tenuation. The experimental attenuation was verified by direct calculation from 
the waveforms. The plane wave attenuation predicted by the Kramers-Kronig rela- 
tion is within 25% the experimentally measured attenuation. The larger than 


expected variations in the plane wave measurements appear to be the result of in- 
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Figure 3-14 Cork-Nitrile Extensional Attenuation versus Frequency 
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strument malfunction. The waveforms were recorded by a digital waveform re- 
corder which appears to have had a slight oscillation in the trigger delay function. 
This would insert a constant phase shift error for a particular center frequency 
group causing the phase speed and attenuation results to be slightly shifted relative 
to other center frequency groups. Review of the extensional phase speed and lon- 
gitudinal phase speed shows the same anistropic behavior as the corprene. Micro- 
scopic inspection of the cork-nitrile shows a preferential alignment of the cork 
particles in the extensional wave propagation direction. 

To summarize our experiences with the measurement methods employed: 1. 
The panel measurements consistently provided excellent result with proper overlap 
between individual center frequency groups (with one exception attributed to 
instrument malfunction). 2. The Young’s modulus method, although appearing 
simple, was in reality a difficult measurement to obtain reasonable data from for 
two reasons: (1) the closed cell rubbers were highly attenuating, which required 
locating the data samples too close together, which introduced large percentage 
errors in the distance measurement between data collection locations. (2) the 
Poisson’s ratios for all of the sample materials was unexpectedly low which reduced 
the amplitude of the displacement recorded by the phonograph needle making the 
signals susceptible to noise. It appears the Kramers-Kronig relation agrees well 
with experimental data for isotropic materials but not so well for slightly 


anisotropic materials. 
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Scattering From the Simulated Infinite Cylinder 


Scattering measurements for the uncoated cylindrical Shell were the first 
measurements performed. The frequency range covered by these measurements 
was 5 kHz to 65 kHz. This experiment served as a check of the experimental 
method since similar measurements are well documented. The scattered signal for 
all ka values was inverted relative to the incident signal. The plot of the 
experimentally determined form function versus ka, where k is the wave number in 
water and a is the outer radius of the shell, is plotted in figure 3-17. The different 
point markers and connecting lines indicate different center frequency groups of 
the incident signal. The theoretical calculation of the form function is plotted as 
the solid line superimposed over the experimental data. The agreement between 
the experimental data and the theory is generally excellent. Agreement is poor in 
the region between ka values of 1 and 2 due to an unexplained resonance present 
in the experimental data at a ka equal to 1.5. There should not be any resonances 
in this ka region hence some undetermined experimental factor must be involved. 
Measurements in this ka region (frequencies of 10 kHz and less) are potentially 
Suspect due to the physical size of the measurement tank which makes the 
Separation of the scattered signal from the multiple background echos 
progressively more difficult at low frequencies. Both curves indicate the presence 


of the first symmetric Lamb mode which causes the first sharp dip in the form 
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function at a ka of 3.7. The second circumnavigation of this wave is observed at a 
ka of 7.45 resulting in a Aka of 3.75. The theoretical calculations assume a Young’s 
modulus of 2.2x10!2 dynes/cmZ2. The first resonance dip occurs when the wave- 
length of the first symmetric Lamb mode is equal to the circumference of the shell. 
The corresponding value of ka that the resonance dip occurs at is determined 


below. 





Yep” | ae = 5.53x10° — (326) 
pl = VY) 


on Be p 
ka = Tou oo a = 3.69 (3-7) 


The theoretical value of the ka location of the first symmetric mode 1s slightly 
lower than the experimental ka value indicating a small error in using the tabulated 
values of either the Young’s modulus or other constants. This slight deviation in 
the material constants most likely resulted from the arn eee processes. 
Since facilities were not available to actually measure the elastics moduli of the 
Stainless steel and the deviation is small, the tabulated values will be used. The 
form function versus ka plot also indicated the first antisymmetric mode will not be 
excited for ka values less than 10. (This observation is consistent with the theory 
since for the thin shell the phase speed of the first antisymmetric Lamb mode is 


significantly less than the propagation speed of the water.) 
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The first coating experiment performed was with the shell coated with 1/16 
inch corprene. Figure 3-18 shows the theoretical form function, calculated by the 
procedure derived in Appendix D, as the solid line and the experimental data as 
points. The theoretical calculation assumes a Poisson’s ratio of 0. The agreement 
of the theory with the experimental data is excellent. The form function graph 
does indicated a slight ka misalignment of the theoretical values with the 
experimental values of the form function similar to that observed for the uncoated 
shell. The theoretical and experimental form functions show resonance dips at 
approximate ka values of 1, 2, 3.7, and 4.2. The experimental data for ka values 
less than 2 is sparce. Two significant observations regarding the effect of the thin 
coating are: 

1. the magnitude of the form function is reduced by the coating for ka values 
greater than 4 but, for ka values less than 4, the form function has an overall 
greater magnitude. 

2. the resonance dip caused by the second circumnavigation of the first 
symmetric mode has been completely eliminated. 

Figure 3-19 compares the theoretical form functions for the 1/16 inch corprene 
coating assuming a poisson’s ratio of .4 and 0. For ka values less than 3.5, there is 
little difference. Substantial difference occurs for ka values greater than 3.5. 
Specifically, there is a larger and shifted resonace at a ka of 4.8 and a second 
resonace, mostlikely the second travel of the first symmetric wave at an 


approximately ka of 7.5. These additional resonances are not substantiated by the 
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experimental data. Due to excellent soreerienl between the theoretical and 
experimental values, Poisson’s ratio will be assumed to be 0. 

In order to investigate the effect of thickness on coating performance, an 
additional .25 inches (.635 cm) of corprene was added to the 1/16 corprene coating 
(total corprene coating of 5/16 inches). Figure 3-20 shows the theoretical form 
function as the solid line and the experimental results as points. The agreement is 
excellent between ka values of 1.5 and 8.5. The experimental data shows a 
decrease of approximately 0.5 dB at ka=9.0. This decrease is probably caused by 
experimental error. In the higher frequency range, the compliant coating is 
becoming more attenuating, therefore decoupling sound from the shell more 
efficiently, and theoretically should be approaching the pressure release boundary 
condition. Expernmental errors that could cause a reduced scattered signal at this 
relatively high frequency (56 kHz) are: 

1. the shell may not be completely vertical. 

2. slight misalignment of the shell, 8103 hydrophone, and F-33 transducer. 

3. possible scattering of the incident signal by the 8103 hydrophone and its 
support tube. 

Significant deterioration of the experimental data occurs below a ka value of 1.5. 
This result is consistent with all previous measurements and is attributed to the fre- 
quency range and measurement facilities. The effect of the additional coating 
material is to shift the resonances to lower ka values and to spread the resonance 


dip over a wider ka range (remove the sharpness of the dip). The sharp resonance 
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dip that occurs at a ka of 3.7 on the thinner corprene coating is shifted to a ka 
value of 2.5. The dip also covers a Aka range of 2 whereas the resonance width of 
the thinner material covers a Aka of .5. The depth of the dip is slightly reduced for 
the thicker coating but not as much as expected. The relative peak that occurs at a 
ka of 4.2 for the thinner coating is shifted to a slightly lower ka value and its 
magnitude is increased from 0.88 to 0.93. Noticeably missing from the thicker 
coating plot is the extensive resonance structure that exist for the thinner coating at 
ka values less than 2. Another unanticipated result is that the magnitude of the 
form function after a ka of 7 1s approximately 1 dB lower than the thinner coating 
form function. The peak that exist at the approximate ka value of 4 has reversed 
relative magnitude positions with the flat portion of form function (for ka values 
greater than 7). If a Poisson’s ratio of .4 is utilized in the theoretical calculations, 
| the resonance dip at ka=2.7 and peak at a ka=4.3 are significantly larger than 
| indicated by the experimental data. The thicker corprene layer has therefore 
_ reinforced the evidence for a small Poisson’s ratio. 

For both corprene thicknesses, the scattered signal was inverted relative to the 
incident signal. 

Figure 3-21 shows the form function plots for various thicknesses of corprene 
coatings. It is particularly apparent that for even large thicknesses of corprene (36 


_ times the shell thickness), the resonance structure is present for ka values less than 


0. 
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The final infinite cylinder coating experiment was to coat the shell with .25 
inches of closed cell nitrile. Figure 3-22 shows the results of the experimental. 
_ values (points) plotted with the theoretical form function (solid line). Again, the 
agreement between the theoretical and experimental values is excellent. The 
_ experimental! data shows some Structure in the 3 to 5 ka range that is not indicated 
in the theoretical form function. Two plausible explanations for this additional 
structure are: 

1. seams in the coating were not sufficiently closed to prevent partial 
excitation of the shell. 

2. when the coating was placed on the shell, the curvature caused the coating 
to be slightly stretched to ensure a smooth fit over the shell. The stretching may 
have caused the air bubble in the coating to elongate in the tangential direction 
and shrink in the radial direction causing an effective reduction in the percentage 
of air, relative to a wavelength, between the fluid and the shell. 

As ka values become smaller, the theoretical form function increases rapidly 
which is characteristic of the perfect pressure release cylinder. Ata ka of 
approximately .2, the form function peaks and approaches 0 as ka approaches 0. 
This behavior is consistent with physical reality since at some small ka (large wave- 
length to coating thickness ratio), the coating will cease to "look" pressure release 
and the form function will be dominated by the characteristics of the shell, ie. 


approach 0 as ka goes to 0. 
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Figure 3-23 displays the theoretical form functions for two thicknesses of the 
closed cell nitrile and the pressure release cylinder as a function of ka. The form 
function for the .25 inch thick sample of nitrile is within 1 dB (the standard accura- 
cy measure for most scattering experiments) of the form function for the pressure 
release cylinder for ka values equal to and greater than 0.5. The resonance 
i structure of the shell is effectively removed (there are no resonance peaks or dips 
_ that have a greater variation than 1 dB from the baseline form function). The 1/16 
inch nitrile does not produce a pressure release condition for ka values less than 4 
as indicated by the large resonace dips at ka values of 3.7 and 0.8. As the ka values 
increase, both thicknesses of nitrile show an increased pressure release behave as 
expected. 

Figure 3-24 is the scattered waveform for the uncoated shell and figure 3-25 is 
the waveform for the nitrile coated shell for center frequencies of 40 kHz. The 
echo from the uncoated shell is dominated by the specular reflection. The only ob- 
vious difference in the waveforms (coated shell versus uncoated shell) is the 
presence of a small extension of the signal at the end of the specular reflection for 
the uncoated shell. 

Figure 3-26 illustrates that a .25 inch nitrile coating on a thin shell, with an 
inner to outer radius ratio of .99, will decrease the form function of the shell by 1.6 
dB (the form function of the coated shell remains about the same magnitude 
regardless of the shell thickness but, the form function of the shell increases for low 


ka values as the thickness decreases). 





95 


JapuljAD aseajay oinssasg PUB IJIN JO. 
SASSIUAIIY | OM JOJ UOIIUN WIO.; JOWBISy¥VY fRoNIIOIY, EZ-E 94NBI4 


| 


8 L 9 G v 2 Cc L 





N 


Sa 72 0 


L 


a yl SI Zi MN 
NOLONN4 WaO5s 


911 





96 


(Aquonbol 4 
19}U9D ZHAOp) WIOJoARM JaNeosyeE [JOYS ULL PIayeoousy p7-€ ond! -] 


“ol (00s) au 
000} GL8 OGL G29 00°S GLE OGC Gol 


000 


(one) i @ hae 
(S}JOA) EpNyduy 


LO 





97 


(Aouanbas,y J9)U9D 
ZHAOp) WIOJIABAA JaNBISHIV [PYG UY L PIALOD IHIIN Se-t aingt.4 


le (09s) Ouul| 
00 OL GL8 OSZ Gc9 00S GLE OSC GCl 000 


B0°O— 


+O°;O—- 
(S}JOA) SPNy|Gduuy 


00'0 


¥0°O 





98 


Japul[Ad 
QSPIJayY JINSsdig puke [JOYS UY], UO BUNBOD [LIN 
JO SOYoUT C79 JO} UONOUN.] WO] JayBOSyoe jeONIIOOYT, 9Z-¢ JIN] 


b> 


8 L 9 G v i C 


eet en 00 te Gf cececesescucs 2 





ra L 80 90 +0 £20 
NOLLONM4 WeOs 


v1 


Y 





99 


Finite Shell with End On Incident Signal 


The following observations are made from the waveforms and normalized 
scattered pressure versus ka plots for the coated and uncoated finite cylindrical 
shell experiments for axial incidence. The finite shell is coated with .25 inches of 
closed cell neoprene rubber for the coated case. Figures 3-27 and 3-28 are 
representative waveforms for the uncoated and coated cases. The center frequency 
for the waveforms in figures 3-27 and 3-28 is 60 kHz. The significant observations 
are: 

1, the uncoated shell’s specular reflection is in phase with the incident signal 
whereas the coated shell’s specular reflection is inverted relative to the incident 
Signal. 

2. The specular reflections for both cases have similar magnitudes. The 
coating does not have a significant effect on the magnitude of the normalized 
Scattered pressure, although the coated shell’s normalized scattered pressure is 
consistently slightly larger than the uncoated normalized pressure, for ka values 
less than 11 as indicated in figure 3-29. The reason for this slight magnitude 
difference is that the coated cylinder is physically larger due to the .25 inch in- 
crease in radius. This 12.5% increase in radius affects the scattered field in a 
quadratic manner whereas the ka scaling only accounts for the increased radius lin- 


early. For ka values greater than 11, the neoprene coating eliminates the severe 
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resonance structure present in the uncoated case. 

3. For the uncoated shell, a second scattered signal arrives at 320 
microseconds after the specular reflection. This time difference corresponds to the 
time necessary to travel the length of the shell at the longitudinal wave speed of 
the stainless steel (68 seconds) plus the time to travel the length of the shell at the 
propagation speed of the water (251 wseconds). The difference between the 
experimentally observer time delay and the calculated time delay is less than 5%. 
The actual path traveled is unknown. In figure 3-30, the center frequency is 45 
kHz, the second scattered signal shape is different from figure 3-28 and it appears 
the signal is actually a combination of signals indicating possible multiple 
reflections are occurring in the end caps. 

4. Figure 3-31 compares the coated shell’s normalized scattered pressure 
versus ka with the normalized scattered pressure determined by the SHIP program 
for the ideal pressure release finite cylinder. Both the coated and uncoated shells’ 
normalized pressures follow the trend of and are within 2.5 dB of the SHIP results 
for ka values less than 11. At the larger ka values, the uncoated shell’s normalized 
pressure significantly departs from the pressure release structure predicted by the 
SHIP program. The SHIP program values were verified by the CHIEF program 
and were found to be within 1% at a ka of 2.0 and within 3% at a ka of 3.0. The 
CHIEF program employed a 480 sub-division scheme. Verification at larger ka 
values was not performed due to the large amount of computer memory 


(proportional to the number of subdivisions squared) and the long computational 
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time required by CHIEF. For the equivalent accuracy (as indicated above), the 
CHIEF program required approximately 5 hours computer time (cpu) versus 
approximately 2 minutes of computer time for the SHIP program on a VAX 750 
computer. 

Although the coated shell’s normalized pressure does show some Structure, 
particularly at the low ka values, it is unclear if this results from acoustic transmis- 
sion through the layer and excitation of the shell or from acoustic excitation of 
shell via the seams of the coating. In either case, the resonance effects are substan- 


tially damped. 


Finite Shell with Normal (Side) Incident Signal 


The following observations are made from the waveforms and normalized 
Scattered pressure versus ka plots for the coated and uncoated finite cylindrical 
Shell experiments for normal incidence. The finite shell is the same shell and 
coating used for the end on or axial incidence case. Figures 3-32 and 3-33 are 
representative waveforms for the uncoated and coated cases. The center frequency 
for the waveforms in figures 3-32 and 3-33 is 40 kHz. The observations are: 

1. The uncoated shell’s specular reflection is in phase with the incident signal 
whereas the coated shell’s specular reflection is inverted relative to the incident - 


signal. 
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2. For the uncoated shell, the two types of circumferential waves are 
observed. The first circumferential wave observed is the first symmetric Lamb 


mode which travels with a phase and group speed of | 


ue | 2 7 (ee) 
p(l-v) 


The difference in arrival times of subsequent circumnavigations of the shell is 


it 
an 
qn 
Gd) 
x< 
fmm 
© 


equal to the circumference of the shell divided by the phase speed (60 yuseconds). 
The first travel of the symmetric wave is contained in the specular reflection and 
can not be directly observed. Subsequent travels of this wave is observable in 
figure 3-34 which has a 20 kHz center frequency. The second circumferential wave 
observed is the first antisymmetric Lamb mode. The phase speed is determined 


from equation (1-8) 


: Tae wh 
va ne ae wn 
7 Spall ay) 


This equation for phase speed neglects any curvature effects. Viktorov?® 
demonstrates a curvature correction that relates the phase speed as calculated in 
the center of the shell (by the above equation) to the phase speed at the outer 
surface by 
h 
> +a 


V = Vo = 1.047 V (3-9) 
p, surface a p p 
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where a is the outer radius of the shell. For the 40kHz center frequency incident 
signal, the curvature corrected phase speed for the first antisymmetric Lamb mode 
is 1.53x10° cm/sec. In order for the wave to be generated in the shell, the phase 
speed of the wave must be greater than the propagation speed of the water. Fora 
30 kHZ center frequency, the first antisymmetric mode is still observed even 
though the calculated phase speed is less than the propagation speed of the water. 
Three reason for this are: 

a. The elastic constants for the stainless steel shell are actually higher than 
the tabulated values, possibly due to the fabrication process. 

b. The sound propagation speed for the actual measurement facilities needs 
to be more accurately determined (actual temperature corrected speed is 1.48x10 
cm/sec). 

c. The standard plate analysis is not sufficiently accurate for this shell 
experiment, since the shell is fluid loaded on one side and air loaded on the other. 
_ In addition, end effects caused by the thick end caps are not accounted for. 

Figure 3-35 is the waveform for the coated shell for an incident signal center 
frequency of 40kHz. Particularly noticeable is the lack of the antisymmetric mode 
signal after the specular reflection. 

3. In figure 3-36, the normalized scattered pressure is plotted against 
increasing ka values for both the coated and uncoated cases. The normalized 


Scattered pressure for the coated shell is approximately three times larger than the 





WV 


(Aouonbas.J 19}U9D Z}y OZ) 
QOUIPIOUT [PULION JOJ BUZIGODNY JO S9YDUT SZ" YUM 


parwod (6° =2/q) []2Y4S aU] YON P Woy jeUdIg palayeog CE-¢ andl 


(SpuUCdas) FWLL 
0'9 O'S O'v oa O'~ Ol 0'0 


QQ0°0 SO'0O—- OPO = 


(S}JOA) SGNLNdNY 


S00 





13 


(1l9YUS poweoous) 
SI OUI] poNod ‘suaidoaN JO SaydU! C70 YIM pajwory 
[24S SI aur] POS) aouaplouy [PUION YIM (6'=B/Q) 1194S 
SNUL YOULL JOJ BY SNSIDA JINSSIIg PIBVPIG Pozi[VULION OE-¢ aIndt.] 


SH 
caesaiadgwniadh o6eb bg 298 G¢ + € Zt O 


2 ° 
i : i Oo 
3 Cad fg 
fi AVS 
23 ‘s . 
3 : e / 
“ : 37% 
Fi peas: 
ty it F3 or O 
02 se : : rat 
s e ‘ : j \' ; 
J ad : 
ry : ry ? 
Du ? i WN 
° s Wt’ ; \ 
a ai : ~ Oo 
Wl f On 


< 
a 
a 


(2°0 


add Pa ({D9s)q 





114 


(AALHO SI oury poyseq ‘oustdoan jo soyoul ¢7'— Yim payvor 

[24S SE OUT] POS) JopurfADr oul asvojay siNssotd 

JO} INdINC AAIHO pue suapiouy [LULON YIM (6° = B/G) [94S 
IUULY YOULL, JOJ BY SNSIIA IINSSIId PIJIIeIG PIZI[PULION LE-¢ d1NB1-] 


aid 


| ~& 
Pi. 
P al 
af 
<a of 
-~—_ 
o 
a -* 
-* 
oF 
~a&” 
ao” 
lad 
o 
-* 
y 126 
wo 
Pe 
ad 
” 
aot 
ae 2 
o 
oo” 3 
-* 
o 
- 
a” 
~°°8 
~~ e | 

-" a é 
= a 


GO'0 


GLO Ol'O 
(J24)q/Q09S)d 


770 





115 


normalized scattered pressure for the uncoated shell for ka values less than 4.0. 
The larger scattered pressure is partially accounted for by the coated shell being 
slightly larger, physically than the uncoated shell due to the addition of the .25 
inches of neoprene even though the ka values are adjusted to account for this. 

4. For ka values greater than 6, the normalized scattered pressure of the 
coated shell shows a significant reduction in Structure. It appears that the 
resonances are averaged to a slowly varying mean value of the uncoated shell’s 
normalized scattered pressure. 

5, The experimentally measured normalized scattered pressure for the coated 
Shell is compared with the ideal pressure release scattered pressure calculated by 
the CHIEF program in figure 3-37 for ka values less than 4.0. All of the data, with 
one exception, is within 1 dB of the calculated pressure. The experimental 
pressure oscillates about the pressure release value indicating the average pressure 
over the ka range approximates the pressure release condition. The CHIEF 
program was run with 480 subdivision and is accurate for a ka value less than 4. 
Subdivision schemes necessary for ka values greater than 4.0 were computer 


memory and time limited. 
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CHAPTER IV 


Conclusions and Recommendations 


The primary objective was to determine if a shell, thick or thin, can be made 
to scatter as a pressure release structure by the application of a viscoelastic coating. 
In order to achieve this boundary condition, the coating must cause the following 
Briects: 

1. The coating must cause the pressure on the surface of the target to be sig- 
nificantly smaller than the incident pressure. The physical results of the signifi- 
cantly smaller surface pressure are that the scattered signal must be inverted 
relative to the incident signal and the magnitude of the scattered pressure must be 
equal to the magnitude of the theoretical pressure determined for the pressure 
release solution. 

2. The shell must be acoustically decoupled from the scattering problem, 1.e. 
the resonance effects must be removed and the coating must not exhibit elastic 
effects of its own. 

The secondary objective was to determine how thick the coating must be in order 


to achieve the primary objective. 
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In order to determine the thickness of coating required, the viscoelastic 
constants of the coating materials had to be determined. For the closed cell 
rubbers the plane wave modulus was determined for the frequency range of 15 kHz 
to 30 kHz and the Young’s modulus for 500 Hz to 3 kHz. The Young’s modulus 
was then extrapolated to higher frequencies in order to estimate the Poisson’s ratio 
and the shear modulus. The Young’s modulus extrapolation was not as accurate as 
had been anticipated based on previous observations for other materials. This 
inaccuracy resulted in a wide spread for the derived Poisson’s ratio, (estimated 
between 0 and .2), causing a correspondingly wide spread for the complex 
transverse wave number ( 1.41 to 1.63 times the complex longitudinal wave 
number). There was no measurable difference in the theoretical form functions 
determined by the two different values of the transverse wave numbers. This in- 
Sensitivity of the form function to the variation in the the transverse wave number 
results from: (1) the phase speed of the closed cell rubber is at least an order of 
magnitude smaller than the phase speed of either the shell or the surrounding 
fluid, therefore, the 13% change does not significantly affect the overall speed 
relationships. (2) the attenuation in the closed cell rubber is very large which 
results in an almost complete decoupling of the shell. An increase in the attenua- 
tion caused by increasing Poisson’s ratio from 0 to 0.2 can not decouple the shell 
further. The experimental and theoretical results are in excellent agreement. The 
cork-rubber composites were determined to be slightly anisotropic. This anisotro- 


py is small enough not to invalidate the isotropic model. The Poisson’s ratio was 
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estimated to be very small (material properties are dominated by the cork which 
has a Poisson’s ratio of approximately 0). For the corprene, a Poisson’s ratio of 0 
resulted in excellent agreement of the theory with the experimental results. 

For the (infinite) thin shell, the pressure release boundary condition can be 
achieved by the addition of the .25 inch nitrile coating. The inversion of the 
specular relative to the incident signal can not be directly attributed to the coating 
since the shell’s form function is in the soft background region for the ka range of 
study. For thinner coatings of the nitrile, the pressure release boundary condition 
can only be achieved at higher ka values. 

The corprene does cause a reduction in the resonance effects, but even at 
_ large thicknesses (3.2 cm), the corprene is not as effective as the .635 cm of nitrile. 
_ The overall conclusion is that corprene is not a viable pressure release coating for 
_ low ka values. 

For the thick (finite) shell, the .25 inch neoprene coating is a fair 
approximation to the pressure release surface. For both orientations of the target, 
the neoprene causes the inversion of the specular reflection and a significant sup- 
pression of the resonances. The normalized scattered pressure for the axial inci- 
dence case is within 2.5 dB of the theoretical pressure release value (mostly 
_ Sreater) as calculated by the SHIP radiation program. For the side incidence case, 
the normalized scattered pressure is within 1 dB of the pressure release value as 


_ Calculated by the CHIEF radiation program for ka valued less than 4.0. 
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For all the shell measurements performed, the application of the coating 
significantly increased the amplitude of the scattered pressure relative to the 
uncoated in the low ka range. The .25 inches of nitnle added to the thin shell 
(inner to outer radius ratio of .97) raised the form function 4.5 dB relative to the 
uncoated shell at a ka of 2.0. The .25 inch thick neoprene coating increased the 
normalized scattered pressure of the finite thick (inner to outer radius ratio of .90) 
shell 8.9 dB relative to the uncoated shell for side incidence at a ka of 4.0. (There 
was not a significant change in the magnitude of the normalized scattered pressure 
of the finite shell for end on incidence indicating that the end caps were scattering 
as a rigid body.) As the shell becomes thinner, however, the closed cell coatings 
will actually lower the for function. 

All of this data was taken over a frequency range of 15 kHz to 70 kHz. It is 
' not expected that these results would scale well to any other frequency range 
(typical ka scaling) because the performance of the coating depends on the size of 
the air bubbles in the material relative to the acoustic wavelength. The size of the 
_ air bubbles in the materials is fixed and therefore, for any other frequency range, 
the ratio of bubble size to wavelength would change. 

The materials in this study are not well suited for deep submergence 
applications. Corprene has been shown to have an increase in its longitudinal wave 
speed of 300% due to pressure effects (see Horsely and Thompson! 16), The 
neoprene closed cell rubber was left under 15 feet of water for several weeks which 


resulted in an unrecoverable reduction in volume of approximately 50%. 
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The acoustics community has long considered corprene as the benchmark 
pressure release material for underwater scattering. This study demonstrates both 
experimentally and analytically that corprene is, at best, a fair pressure release 
material and Is substantially outperformed by the closed cell rubbers. The 
tendency of the closed cell rubber to degrade with submergence time (the air 
leaches out of the rubber) can be prevented by potting the closed cell rubber ina 
material acoustically matched with the water (i.e rho-c rubber). It is also apparent 
from this study that the thickness of the shell is an important variable when deter- 

_ mining the thickness of the coating required to achieve the pressure release 
condition (.25 inches of the closed cell rubber decouples the thin shell more effi- 
ciently than the thick shell). 

The following continuing work is recommended: 

1. The development of techniques which would allow the shear modulus to be 
measured in the frequency range of interest. 

_ 2. All theoretical models of the materials which could predict the shear and bulk 
moduli failed to predict accurate values. A model that allows air concentrations of 
approximately 80 volume percent is necessary. 

3. A numerical computer code that allows the evaluation of the scattering problem 
for a finite structure that takes into account the elasticity of the viscoelastic 


coating, the shell and the end caps. 
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APPENDIX A 


Derivation of the Scattered Field and Surface Velocity From an Infinite Pressure 
Release Cylinder for an Incident Plane Wave 


In this appendix the scattered pressure from an infinite, pressure release 
right circular cylinder is determined for a normally incident plane wave. The 
surface boundary condition (pressure release) defines the acoustic pressure on the 
surface to be zero; this allows the actual material properties of the cylinder to be 
neglected. The derivation is general in that the incident plane wave is taken to 
approach the cylinder at an arbitary oblique angle, a, to the center axis of the cylin- 
_ der (see figure A-1). The solution will be reduced to the special case of interest i.e. 
normal incidence by setting @ equal to zero. 

The governing equation for the region exterior to the cylinder (r>a) is the 
Scalar wave equation 


aP(r,t) 


até 


vp(r,t) = ~ (A-1) 
Cc 


A harmonic time dependence, exp(-twt), is assumed which reduces the wave equa- 


tion to the scalar Helmholtz equation 
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Poexp(LK-X) 





Figure A-1 Geometry for Infinite Pressure Release Cylinder 
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Lae) 


yoP(r) + k°P(r) = © where k° = = (A-2) 


C 


In the previous two equations, P(r,t) is the total acoustic pressure and consists of 
the incident pressure, P(r,t)jn¢, and scattered pressure, P(r,t),. The acoustic wave 
speed of the exterior fluid medium (water) is denoted by c. In the Helmholtz 
equation, the time dependence is suppressed. In order to determine the form of 
the scattered pressure, the Helmholtz equation for the geometry of interest (right 
_ circular cylinder) must be investigated. For the infinite right circular cylinder, the 


Helmholtz equation is separable. The pressure is assumed to be of the form 


P(r,@,z) = R(r) 6(@) Z(z) (A-3) 


where r,@,and z are the standard cylindrical coordinates. The Helmholtz equation 


becomes 








a =a0 (A-4) 


After applying the separation of variables technique, the following three equations 


result 


2 
eal bv tl oe ae (A-5a) 
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2 
cs Oy Se (A-5b) 
dé 
2 
ene » K 7 = 0 (A-5c) 
2 Z 
dz 
2 a 2 2 W W 
kK = — —e= k_ + k_ where k =-—cos(a) and k_ = — sin(a)] 
<A r Z ig ame AINE 


The general solution for the radial equation, R(r), is a combination of Bessel 
functions of the first and second kind (or Hankel functions of the first and second 
kind). The solution for 6(@) iscos(mé@) and sin(mé), where m must be an integer 
for the solution to be continuous at @ = 0 = 2m. Due to the symmetry, the 8 
function must be even with respect to the @ = O which eliminates the sin(m@) form. 
The general solution for Z is a combination of exp(-ck,z) and exp(tk,z). The 
second solution of Z is eliminated by requiring the solution to be finite as z 
approaches infinity. 


For the scattering problem, the incident wave can be written 


Pinc(r,9,Z) = Po exp(-ck-r) = Pg exp(-vkz)) om Em cos(m@) Jm(kyr) 
m=0 


(A-6) 
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where €,, is the Neuman factor (€,, = 1 form = 0 and €,, = 2 form > Q), 


Jm(k;r) is the cylindrical Bessel function of the first kind, and k is the wave 


propagation vector. 


The total pressure is the sum of the incident pressure and the scattered 


pressure. The scattered pressure is assumed to be of the form 


© 


Po(r,8,z) = Po exp(-tkz2)) uM €, cos(m@) Am Hia(kyr) (A-7) 


m=0 


In order to determine the constant Ap, the pressure release surface bounda- 


ry condition 1s employed. 


© 


Ptotay(a,9,z) = Po exp(-tkzz)) a Em cos(mé@) (An H1(kpa) + Jmn(Kya) ) 


m=0 
= 0 (A-8) 


The orthogonality of the cosine function is utilized by multiplying equation (A-8) 
by cos(n@) and integrating with respect to @ for @ between 0 and 21. 


@ 


on 
Po exp(-ekzz)) eM €n | AD H1 (kya) + Jn(Kra) | | cos(m@) cos(né@) dé = 0 


0 
m=0 


(A-9) 
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After simplification, the constant A,, is determined to be 


J (ka) 
a (A-10) 


H (k_a) 


r 


The total pressure is determined by substituting equation (A-10) into 


equation (A-7). 


i: 0] 


Peotay(r,9,z) = Po exp(-vkzz) uM €, cx] Um(kpa) - 


m=0 


m Yr 
hk) (A-11) 
m 


The actual quantities of interest are the radial surface velocity and the far- 
field scattered pressure for the normally incident (a2 = 0) plane wave. The normal 
incidence condition removes the z dependence from the solution since k, = 0 and 


k, becomes equal to k = w/c. The radial component of the velocity follows from 


Euler’s equation 


VP(r,@,t) = = ae li on 


ar “vr 36°90 ~  ? at - -p(-w)| Vie + Wye, | 





the radial velocity becomes V_ = - 





127 


Substituting equation (A-11) for the total pressure, P, then se becomes 











-Po d 
Vy(a,8) = pm+] €m cos(mé) | Jm(ka)  - 
pc d(ka) 
m=0 
Jm(ka)  d 
a H1.(ka) | (A-13) 
Hl (ka) d(ka) 


The surface velocity, equation (A-13), is simplified by using the Wronskian formula 


d 


] d ] 
W(J,H) = H (ka) d(ka) Jb ka) = J tka) d(ka) H 


m 





f-<2 
(ka) = - m ka (A-14) 
the surface velocity becomes 
-Po 2 ] 


Cue mes mG) 
pc naka H1 (ka) 








(A-15) 


m=0 


The scattered pressure is derived by subtracting the incident pressure from 


the total pressure yielding 
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© 


Jm(ka) 





P.(a,@) = Po uM ¢€,, a] - ou) (A-16) 


HI (ka) 
m=0 

The general form for the nondimensional, normalized far field pressure is 

P(r, 4) 


| 3 ukr 
ols, © 


The far-field scattered pressure is obtained by taking the limit of the Hankel 


= Function(ka,@) (A-17) 


function as the argument asymptotically approaches infinity in equation (A-16). In 
Arfken!17 the asymptotic form of the Hankel function of the first kind is shown to 


be 


1 2 UE = (m+ 5) 5) 
ME - | a (P(g) + 4 OE) ver 


for -m<arg(€)<2n (for this particular derivation € is real and arg(€)=0) 


(A-18) 








iz? 


In the limit as € > o, P(g) + 1 and Q(¢) + 0 resulting in 


1 ° be -i(m + a > 
He) = | a = z (A-19) 


The normalized far-field scattered pressure becomes 





P.(r,@) Jm(ka) 2 1) 1% 
—_——_——————-__ = - uM €, cos(mé) — exp/-ejm+—-] - 
H! (ka) tka 2 uid 
a 
p | — e-tkr m=0 (A-20) 
0 r 


The far field cylindrical form function is defined to be the square root of 2 times 


the magnitude of equation (A-20). 





(A-21) 








The right side of the equation (A-20) is numerically calculated. The Bessel 
functions are calculated using the recurrence and normalization methods described 
in Abramowitz and Stegun!18 and Hitchcock!19. Satisfactory verification of the 
Bessel functions with published tables for test arguments was obtained. (The 


method is also valid for Bessel functions of complex arguments which are required 
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in Appendix C.) The infinite series is terminated when the ratio of the next term in 


the series to the total sum is less than 10-6, 


Determination of the Low Frequency Limit 


Only the first term of the series of equation (A-20) needs to be retained in 
order to determine the low frequency limit of the normalized scattered pressure. 


The normalized backscattered (6 = ™) pressure is approximated as 


n 
P(r, 8) J (ka) l ) 


= Limit ad ——- 23 (A-22) 
a tkr 
Po | e 


] nm Ka 
Hy (ka) 
where the Limit, Ja (ka) = ] and where the 


i . 2 ka ; 
Limit, No(ka) = —— { In i] + vf (y is Euler’s constant) 


The far-field backscattered cylindrical form function becomes 


a | See (A-23) 


te ST) ais 
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Determination of the High Frequency Limit or Specular Reflection 


The following specular reflection derivation for the cylinder parallels the der- 
ivation for the sphere in Urick120. There are four basic assumptions required for 
this method: 

1, the scatterer is perfect in shape, 

2. the magnitude of the reflection coefficient equal to unity, 

3. the scatterer is immoveable, 

4. the radii of curvature of the scatterer is larger compared to 

the acoustic wavelength (ka > >1 ). 

The approach of this method is to relate the energy contained in a small acoustic 
beam which 1s incident on a small area, dS, of the curved surface to the spreading 
of that incident energy upon reflection at a unit distance. For this derivation, the 
cylinder is considered to be infinite with the center axis normal to the wave 
propagation vector. The radius of the cylinder is a. 

As shown in Unick, the acoustic center, Q, is located a distance of a/2 from 
the geometric center, C. If the waves were to penetrate the surface, they would be 
focused at the acoustic center, Q. Upon reflection from the surface area, the 
reflected waves appear to originate from a line source located at the acoustic 
center, Q. The specular reflection far-field form function reduces to the constant 


value of 1.0. 
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P(r, 8) j 
f= —+j{—.- | s wie = 1.0 (A-26) 


The normalized far-field backscattered form function as a function of ka is 
plotted in figure (A-2). Figure (A-2) shows the backscattered pressure rapidly 
approaches the geometrical limit and is within 5 percent of that limit by a ka of ap- 
proximately 6. For the low frequency limit, the backscattered pressure is within 1 


percent of the low frequency limit at a ka of approximatedly .09. 
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APPENDIX B 


Pressure Release Sphere 


The solution for the scattered far-field pressure of the ideal pressure release 
sphere follows the standard approach of expressing the field pressure as an infinite 
series of spherical harmonics and solving for the series coefficients. The geometry 
of the problem 1s illustrated in figure B-1 where the incoming plane wave is trav- 
eling in the + Z direction. Standard spherical coordinates are used. Due to the 
symmetry of the problem, the solution is independent of ¢. The exp(-twt) time de- 
pendence is suppressed. 

The spherical harmonic expression for the incident plane wave traveling along 
the polar axis is 

w 
Pinc(9,%) = Po ) ten + 1) eM P(cosé@) jn(kr) (B-1) 


m=0 


where P,,,(cos@) is the Legendre polynomial of order m and j,,(kr) is the spherical 
Bessel function of the first kind of order m. 
The total pressure, po, is the sum of the incident pressure, Pinc, and the 


Scattered pressure, Pccat- 
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Poexp(LKX) 


Sphere of radius a 


Figure B-1 Geometry for Pressure Release Sphere 
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Prot (0298) a Pec lt9s9) + p estoy) (32) 


scat 
The scattered pressure can be expressed as the following general spherical 


harmonic infinite series 


~ m 
Pocat(r, 4,9) -) ) YOm,n(9,¢) G°m n(kr) where o is +1 or -1 (8-3) 
m=0 n=0 


and Y°rm n(9,¢) is the spherical harmonic and q°;y (kr) is a combination of jp)(kr) 
and h?,,,(kr) which are determined from the boundary conditions. Since the 
problem is independent of ¢, Y°rm (9,¢) reduces to P_,(cos@). The pressure 
release condition causes the incident wave to be completely reflected resulting in 


only an outward traveling scattered wave. 


Gn n (KP) : h_ (kr) (B-4) 


With simplification, the scattered pressure becomes 


@ 


pscat(r,@) “) hn Pr(cosé) hi, (kr) (B-5) 


m=0 





oy 


The total surface pressure is zero and with the substitution of equations (B-1) and 
(B-5), equation (B-2) becomes 
0 = | Po(2m + 1) o™ Pr(cosé@) jm(ka) + Am Pm(cosé) h1 (ka) (B-6) 


=0 


In order to determine the constants Aj, the orthogonality property of the 
Legendre polynomial 1s utilized by multiplying equation (B-6) by P,(cosé@) and 


integrating with respect to d(cos@) from -1 to 1. 


@ 
] 
0 = Po(2m + 1)e™j,(ka) + a) | Pm(cosé) Pyr(cosé) d(cosé@) 
-] 
m=0 
: (B-7) 
Simplification yields the series constant 
m , 
Pt (2m + 1) J,4ka) 
Bs ac (B-8) 
h (ka) 
m 


Substitution of equation (B-8) into expression for the scattered pressure, equation 


(B-5), results in 
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= Po ym (2m + 1) Jm(ka) 
Beeat(r,@) = ———_—_—__—__—_———  P,,(cos@) hl, (kr) (B-9) 
hl (ka) 


The far-field scattered pressure is obtained by utilizing the asymptotic form of the 


spherical Hankel function for the argument approaching infinity 


LZ 
e 


Z 





ae ] m+] 
Vimit h(Z) = (-t) (B-10) 


In addition, the far-field scattered pressure should be of the standard form 


Pee(t8) = py =] f(ka,0) es" (B-11) 


where f(ka,@) is the far-field spherical form function. The resulting normalized far- 


field scattered pressure is 


Pscat(r, 4) (2m + 1 ) jm(ka) 
Se pcos), (B12) 


a ka hi (ka) 
p ~ exp(ckr) 
0 


r m=0 
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APPENDIX C 


Modification of the Naval Ocean Systems Center (NOSC) CHIEF Program 


With NOSC Technical Document 970/74, Benthein and Barach released a 
modern fortran program to compute the acoustic radiation from an arbitrary- 
shaped body using the CHIEF method developed by Schenck/!, The radiation 
| program computes both the near-field and far-field pressures of a closed surface, S, 
provided the normal surface velocity distribution is known. In particular, the 
program first calculates the unknown surface pressure by solving the surface 


Helmholtz integral equation 


EKY UX, G) 


9 | a tkr(x,0) | A 
2nP(x) = P(a) an T(x,0) + LWP V(o) ~ ¥(x,0) _ dS(o) (C-1) 
S 


The integral is solved numerically by subdividing the surface into N, subregions 
and assuming that the surface pressures and normal velocities are constant in each 
Subregion. The result of this approximation is a system of linear algebraic 


equations. The algebraic equations are generated from the following 
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Ns 
exp(-tkr(&,9) 
2nPa(x) - — a 7 | dS(o) = 
YiSmI9 
n=] Sp 
Ng 
exp(-ckr(&p,0) 
Lwp Vn [ee | dS(o) (G22) 
r(€,,0) 
n=] Sn 


where €,, is some reference point on S,,. In matrix notation the previous equation 
can be written as 


Ap=Bvy (C-3) 


where 


| eo kr (fo) | 
Bn = lwp “rE o) dS(a) (C-4) 
S 


and p = (P,,) and v = (V,,) are vector of the surface pressure and surface 


velocities respectively. The surface integrals is A and B are solved by Gaussian 
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quadrature. Matrix equation (C-3) does not yield a unique solution at certain 
characteristic frequencies. The CHIEF method overcomes this uniqueness 
problem by forcing additional constraints on the surface pressure. In particular, 
the interior Helmholtz equation must also be satisfied by the surface pressures and 
velocities for interior locations at all frequencies. This program allows interior 


points of S to be input which are related by 


: a tkr(x,0 ) | ot kr(x,9) 
P(o) pu ¥(x,0) + LWP V(o) “r(x,0) dS (a) (C=3) 
0 
S 
mM = No ae We + N. where N. is the number of interior points. 


Equation (C-3) is still valid but with the coefficients modified to 


-ukr(§ 0) 
nb - | | ah | eee | dS(o) m=l1toN 
mn dn ries. o> S 
S 0 m 


A = 
‘ih 3 oe bkr (x 0) 
| | an | a ala dS(o) m=N toN+N, 
n i Xoo) S eee 
0 m 
S 
| oe tkr(€_ a) | 
LWP —r(é,0) dS(o) m= 1 to Ne 
S m 
n 
B= 
mn 
-ukr(x_,o) 
LWP ee | dS(o) m=N toNdN, 
r(x 0) S Son 
S 
n 


(C-6) 
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The matrix equation is now overdetermined and solved by the Householder reduc- 
tion technique. Once the surface pressure is determined, the exterior Helmholtz 


equation is used to determine the field pressure. 





Nn 

] 0 exp(-ckr(x,o)) 
= — Pr — ed dS(o) 
4n dng mee) 
n = ] Sn 
Nn 
LWwp exp(-ukr(x,c)) 
+ Vn | —_—__—_—_———_ | dS(c) (C-7) 

4n r(x,o) 

n=] Sp 


The previous discussion is a summary of the methodology used in the program 
developed by Benthein and Barach with considerable more detail contained in 
NOSC TD 97074. 

The objective of the appendix is to modify this program to determine the 
scattered pressure field from a plane wave incident (Pj,,) on a pressure release 
surface. Unlike the radiation program, the surface pressure is initially defined to 
be zero with the normal surface velocity to be determined. 


The scattering version of the Helmholtz integral is still applicable. 
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r(x,o) r(x,o) 


3 eo tkr(x,0) etkr(x,0) 
4nP(x) = | P(o) an | Too Oo —————— 9 
0 
S 


+ 4n Pax) (C-8) 


Since P(x) is defined to be zero on the surface, the equation (C-8) reduces to 





1 ol kr(x,0) 
PAX) = - a | twp V(o) —r(x,0) dS(o) (C-9) 
S 


The surface geometry and surface subdivision sections of the the CHIEF program 
are utilized to define the surface. The incident pressure 1s calculated at the center 
of each subdivision by exp(-tkX3) assuming the plane wave is traveling in the + X3 
direction (this program assumes an exp(twt) time dependence). Matrix equation 
(C-3) reduces to 


Bv = - P. (C-10) 


from which the normal surface velocities are determined by a standard equation 
solver subroutine. The difficulty of characteristic frequencies still exists but matrix 
equation (C-10) may be overdetermined in a similar manner to the radiation 
problem. Once the normal surface velocities are known, the field pressure 
portions of the CHIEF program are utilized, with the surface pressures set to zero, 


to compute the scattered field. The scattered far field pressure is computed from 
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Ns 
P(x) 
Pee(x) So a d\n Be¢(x,n) | (C-11) 
uk exp(vckR) 
(a n=] 


4x R 


where Bef is the far-field version of B. 

Verification of the program modification is only expeditiously possible for two 
closed pressure release surfaces, i.e. the sphere and the finite cylinder. An exact 
analytical solution is available for scattering from the pressure release sphere and 
is contained in Appendix B. The finite pressure release cylinder with an axially 
incident plane wave can be compared to the proven solution generated by SHIP82. 

The scattering geometry for the pressure release sphere is illustrated in figure 
(C-1). The surface is divided into a total of 144 subdivisions with 12 equiangular 
divisions in both the theta and phi coordinates. If the usual rule-of-thumb of 10 
subdivisions per wavelength is used, this subdivision scheme is anticipated to be ac- 
curate up to a ka of 2.5. Due to symmetny,, the field is only dependent on the r and 


theta coordinates. The analytical solution for the far-field pressure is 





P(r,@) 
Ng 
Po -L(2m+1) jm(ka) 
———______ = )- ——_—_—_—_—_———— P,,(cos(@)) (G12) 
Lk exp(-ckR) hi .(ka) 
——_—_—_____ m=0 


4n R 


where P,,,(cos(@)) is the Legendre polynomial of order m. 
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A comparison between the program and analytical solutions of the 


normalized far-field pressures for ka between 1.0 and 4.0 is tabulated below. 


Comparison of CHIEF Results for Pressure Release Sphere 


-.547 
-.205 
oo 
ek, 
220 
-.295 
PL 
-.247 


Program Solution 
Real 


Imaginary 


581 
IO 
Soy 
57 
ea 
20 
Ss) 
589 


-8.69E-2 


-.535 
-.424 


6.93E-2 


478 
428 


-4.56E-2 


-.534 


Program Magnitude 


Table C-1 


Analytical Solution 
Real 


-.944 
-.204 
Oe 
28 
210 
-.291 
-.515 
-.256 


80 
il 
BS) | 
eo 
Baal 
518 
14 
ro 


Imaginary 
-8.72E-2 


-.533 
-.422 


6.95E-2 


475 
429 


-1.69E-2 


-.443 


Analytical Solution Magnitude 


Error(%) 


0.20 
0.35 
0.37 
OVS 
0.96 
U9 
2.14 
15.04 
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The pressure release program is accurate to within 1 percent over the ka range of 
1.0 to 4.0. The results are somewhat surprising in that accurate results are 
obtained up to a ka value of 4.0 (less than 10% error) which results in approxi- 
mately six subdivision per wavelength. There are no indications that the two 
solutions would not match at larger ka values provided the number of subdivisions 
was appropriately increased. 

The cylinder verification 1s for the plane wave traveling in the -Z direction and 
incident on the flat ends of the cylinder. The subdivision scheme was to divide the 
top into three circumferential rings of equal radial dimension. Each ring was also 
divided into six equiangular sections in the theta direction. The bottom end was 
partitioned in a similar manner. The side of the cylinder was divided into six equal 
Sections in the axial direction and six equiangular sections in the theta direction for 
a total 72 side subdivisions (see figure (C-2)). For the SHIP solution, the top and 
bottom ends of the cylinder were divided into 3 circumferential rings of equal 
radial dimension. The side of the cylinder was divided into six equal 
circumferential rings in the axial direction. Due to the symmetry of the problem, 
the subdivision schemes are equivalent. As with the pressure release sphere, the 
anticipated accuracy is up to a ka of 2.5. The results are tabulated for normalized 


surface velocity in each of the rings for ka=0.5, 1.0 and 2.0. 
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3 


JapulfA) aNULy aseajay 
IINSSII JOJ QUDYIS UOISIAIPGNS puev A1jJaWOIH BulsayeIS 7Z-d sINBI] 


MAaIA GNA 


(QN3 MOVE SALVOIGNI SISSHLNIYVd) 


(ZL GNVG) © GNVG 
(LL GNVS) Z GNVd 


(OL GNVS) | GNVE 


= SAVM LNACIONI 














Oo OonN NN SP WN 





Oo OTN NP WN 


-1.207 
-1.254 
-1.806 
-2.226 
-0.9588 
-0.6741 
-0.4270 
-0.1997 
0.2568 
0.5964 
0.5193 
0.5199 


0.4213 
0.4854 
Ue, 
Lee's) 
129 
1.248 
eS 
e296 
2.049 
1.046 
0.5266 
0.4618 


Table C-2 


Pressure Release Program 


Real Velocity (x10-6) 


Imaginary Velocity (x10- 


Comparison of CHIEF Program With SHIP Program | 


SHIP 


Real Velocity (x10-6) 


-1.203 
-1.256 
-1.808 
-2.219 
-0.9593 
-0.6740 
-0.4269 
-0.1994 
0.2554 
0.5968 
Cian97 
0.5178 


Imaginary Veloci 


0.4207 
0.4869 
0.9953 
1.991 
1.280 
1.249 
1.254 
1.296 
2.046 
1.050 
0.5281 
0.4610 
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Real Veloci 


Pressure Release Program 


-1.033 
-1.081 
-1.595 
-2.040 
-0.8935 
-0.6263 
-0.3906 
-0.1750 
0.2171 
0.4387 
0.3467 
0.3350 


Imaginary Veloci 


-0.2564 
-0.2099 
0.0132 
0.7538 
0.6674 
0.7140 
0.7416 
057525 
1.036 
0.3033 
0.0171 
-0.0303 


x10-6)Real 


HIP 


Velocity (x10- 


-1.030 
-1.082 
-1.597 
-2.034 
-0.8941 
-0.6263 
-0.3905 
-0.1758 
0.2159 
0.4392 
0.3472 
0.3338 


Imaginary Veloc} 


-0.2545 
-0.2092 
0.0151 
0.7537 
0.6676 
0.7141 
0.7516 
0.7525 
1.035 
0.3049 
0.0176 
-0.0297 
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Pressure Release Program 


Real Veloci 


5.538 
-6.100 
-10.606 
-16.442 
-7.109 
-5.215 
<208 
-0.534 
3.075 
2.126 
0.318 
-0.178 


Imaginary Velocity (x10-7) 


-9.867 
-9.346 
-9.135 
Pa Ue 
2.853 
4.439 
5215 
5.230 
5.587 
-0.1290 
-1.395 
-1.544 


SHIP 


Real Velocity (x10-7) 


=e 20 
-6.115 
-10.633 
-16.402 
-7.714 
-5.216 
-2.779 
-0.535 
3.063 
746) ie 
0.319 
-0.175 


Imaginary Veloci 


-9.823 
-9.348 
-9.127 
-2.285 
2.851 
4.439 
3212 
Dice! 
5.583 
-0.1216 
1395 
-1.541 
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APPENDIX D 


ttering From th Infini linder 


In this appendix, the solution for the far-field scattered pressure from an 
infinite, hollow elastic cylinder coated with a viscoelastic solid is derived. The 
elastic cylinder is considered lossless whereas the losses in the viscoelastic solid will 
be accounted for by employing a complex wave number. The method of approach 
is to formulate the displacement potentials for the hollow, infinite elastic cylinder 
for an arbitrary incident plane wave and reduce this solution for the normal inci- 
dence case. This result will then be extended to the layer problem. This solution is 
applicable to all layer and wall thicknesses since no shell approximations are 


utilized. 


Scattering From an Infinite Elastic Cylinder With Arbitrary Incidence 


The geometry for the elastic cylinder is shown in figure D-1 where y is the 
angle of incidence referenced to the Z axis. In the elastic solid, the displacement 


equation of motion for no body forces is 
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Poexp(uK-X) 





Figure D-1 Scattering Grometry for Infinite Cylindrical Shell 
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rv vu + (A + 2n)V(V-U) = p — (D-1) 


where U is the displacement, t is time, p is the material density, and up and dX are 
Lame constants (p is the shear modulus and + 2y is the plane wave or longitudi- 
nal modulus). Utilizing Helmholtz’s theorem, the displacement is written as the 


gradient of a scalar potential function and the curl of a vector potential function. 


U(r,6@,z) = ue + ve, + We a YW+tVxy¥ (D-2) 


Associated with this expression for the displacement, the gauge invariance of the 
transformation of equation (D-2) allows the divergence of the vector potential to 
be an arbitrary function (see Armenakas!21), 


V- ¥ = F(r,t) (D-3) 


Substitution of equation (D-2) into equation (D-1) yields 


2 
wV(VS + Vx H) + (A+ MH) VV» (HV ¥) = pA (WEFT XY) 
dt 


This equation is simplified by using the following vector identities 


V - Vo = v4 and V-Vxw¥=0 


The equation of motion is then reduced to 


2 2 
V/s nv - 9 56 | sox [ev¥-p—z¥]-0 (D-4) 


at ae 
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From equation (D-4) it can be seen that the potential expression for the 
displacement, (D-3), satisfies the equation of motion if the scalar potential, ¢, 
satisfies the scalar wave equation and the vector potential, ¥, satisfies the vector 


wave equation. 


v¢ = agg where ¢? = mame (D-5) 
2 2 d p 
6 dt 
d 

2 ] a°y 2 

V¥= roa where C. oe (D-6) 
CY at 3 


In equation (D-5), Cg is the dilatation or longitudinal wave speed and C, is the 
transverse or Shear wave speed. Each wave equation will be solved using the 
standard separation of variables technique. 


The scalar wave equation, (D-5), in cylindrical coordinates 1s 


9° l1 a ] 9° 9° : 
meee ai) pe ary ¢d (D-7) 
or r or r 06 02 C 


wt ; ; 
(a harmonic time dependence is 


where $(r,8,z,t) = R(r)@(8)Z(z)e 
assumed). The functional form of ¢ is substituted into equation (D-7) 


and after simplification, the equation can be written as 


Sag Se Ee eee (D-8) 





156 


The Z(z) dependence is separated out with ¢¢ being the separation con- 


stant yielding 


MewE 7 = 0 which results in Z(z) = Ae“? + Be 7 | 


(D-9) 


Equation (D-8) is modified with the second separation constant, n, to 


The theta equation separates out to 6" + n°6 = 0 


ea) 


where 6(6) = D cos(@n) + E sin(@n). Since 6(@) must be single valued, 


n must be an integer, n, equal to 1,2,3,... 


The radial equation can be written as 





The most general solution is: 


(p12) 


(D-13) 
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R(r) = GJ far) + HY (ar) : 


The general form for ¢(r,6,z,t) is 


© 


J (Anexp(-eéz) + Brexp(-t€z)(Dpycos(n@) + E_sin(n@)(GpJp_(ar) + HnYp (ar) ) 


n=0 
° exp(-twt) = 0 (D-14) 
The incident pressure wave for an angle of incidence of 7 is 


ot (KeR - wt) 


Berir,0.z,t) = P. 


pa ) u" €, cos(né@) an(— cos(yr)) exp(-ewt) (D-15) 


n=0 


(C is the propagation speed in the fluid) 

For this derivation, ¢ represents the dilatation wave which results from the incident 
wave and, should have a form similar to its exciting force. In particular, the Z de- 
pendence should be the same as the incident wave, i.e. a traveling wave in the Z 
direction, which requires B, = 0. The potential ¢ must be symmetric with respect 


to @ = 0 (even funtion in @) requiring E, = 0. The result is 
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© 


g(r,8,z,t) = exp(r(€z - ut) cos(n@)(GpjJn(ar) + HAYn(ar)) (D-16) 


n=0 


The vector wave equation in component form with an e7!“! time dependence 


iS 








2 * 2 ¥ 2 . 2 Er 
0 a Cant Vata —_ + — 5 e, 
r r 0g Yr r 06 

2 “i 

+Vv e = ¥e+¥, e +¥_ e (D-17) 
zZ 2 ce rr 6 @ ZZ 
S 
The following solutions are assumed: 
¥(r,8,z) = h (r) sin(nd) ace (D-18) 
¥a(r,8,2Z) = hy (r) cos(né) el 62 (D-19) 
LEZ 

¥_(7,8,2) = h(r) Sin(n@) e (D-20) 


The @ functions must be singled valued which requires n to be an integer. The se- 
lection of the sine and cosine term in ¥,, ¥g and ¥, results from the required nature 
of the @ dependence for the longitudinal, torsional, and flexural modes (see 
Graff!22 p.466). The form of coupling between ¥, and Wg in the e, and eg 
components indicates that a sine dependence for ¥, is consistent with a cosine de- 


pendence for ¥g. Additional justification for the trigonometric selection for ¥, will 
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be presented later. The Z vector component results in the scalar wave equation 


for ¥, 


The solution to equation (D-21) is similar to that for ¢(1,@,z) and is 


¥(r,6,2) = ef ($2) sin(nd)(C,J_ (Br) + CY (Br)) 





where B = 7 E 


The two coupled equations for the h,(r) and hg(r) dependences are 


dh an 2 2 
—- + - — + — -nh + 2nh. - h - €h + 
2 2 r 8 Yr r 

dr , oF r 

d*h 1 dh 
i, - n he + énh - hy - € Ne + 
dr r dr r 


(D-21) 
(D-22) 
(D-23) 
a 
—h =0 
c2 r 
S 
(D-24) 
oe 
—h, =0 
me 8 
S 
(D-25) 


Equations (D-24) and (D-25) are then solved simultaneously for 


he and h, resulting in two coupled Bessel’s equations. Subtracting 


(D-25) from (D-24) yields 
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er oe hy =.) SO (D-26) 
2 2 r 8 
dr (lle r 
This Bessel’s equation solution is 
h - h = 20, J J fF") + 2C , Y le (D-27) 
Adding equations (D-24) and (D-25) yields 
2 2 
ee ge MM th th) = 0 (D-28) 
2 2 r 6 
dr Gar r 
The solution to this Bessel’s equation is 
he + h, = 2C. J J (Br) + 2C, ” mnie) (D-29) 


adding and subtracting equations (D-28) and (D-29) gives he and Ne: 


= 
l 


_ CJ, 4Pr) + C, \ ne (Pr) + C 5 -4Pr) + C. Y ,.48r) (D-30a) 


Ne = CJ, {Pr) +C 4’ na hr) - C 5 ne br) - C. Y 4 Pr) (D-30b) 


As shown by Armenakas, et al.121 the gauge invariance yields¥ - ¥ = F(r) 

where F(r) is arbitrary. Meeker and Metzler!23 demonstrate that F(r) can be 
adjusted so as to eliminate one potential function (since there are eight constants 
and six boundary conditions); this is equivalent to setting C5 and C¢ equal to zero. 


The result is 
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fe = -h where h = C494 48r) + CAY ,{er)- (D-24) 


r g r 


The displacement components of equation (D-2) can now be written as 


‘ee | R’ + — h, + éh eoeiijee aun”? (D-25a) 
n U(€z-wt) 

V = as R + ch. - ns | sin(n@) e (D-25b) 

ven -gR - ih’ ant) h FRU OO eae 


; n+] 7 
siny = 0) or he + : he = 0 


For normal incidence, w = 0 (€ 


Substitution of equation (D-24) into the above expression yields 


n+] : n+] _ 
% [ 9% {er 7 sae " 4 | eee une se ‘ ng”? = 


This expression must be valid for all values of n and fr. If n = -1 the above expres- 


sion reduces to 


C,9, (Br) = -C,Y¥, (Br) 


In the limit as Br + 0,C4 must equal 0, and in the limit as Br + ©,C3 must equal 0. 


_ This is only possible is both C3 and Cq are zero; therefore, h, and hg mustalso be 


zero for the normal incidence case. This result is also intuitive since the only con- 


tribution to the Z displacement component, which is zero, is from the curl of h, 


and hg. The gauge invariance allows one potential to be set to zero, therefore the 
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other must also be zero or constant. 

For the normal incidence case, ¥, represents the shear wave in the cylinder 
which must vanish at @ equal to” and 0, and be a maximum at @ equal to 1/2 and 
3n/2. Since the shear wave is antisymmetric with respect to @ equal to 0, the sine 


trigonometric function is the proper selection for the @ function for ¥5. 


Scattering From the Two Layer Cylinder With Normal Incidence 


The normal incidence layer problem is similar to the elastic cylinder (normal 
incidence) derived previously in that two potential functions can fully represent the 
wave behavior in layers. The outer layer is assumed to be viscoelastic with the loss 
factor being represented by a complex wave number. The inner layer is assumed 
to be perfectly elastic which is valid for most metals below the megahertz frequen- 
cy region. The method of solution is to represent the stresses and displacements in 
the four regions and solve the boundary value problem. The geometry is 
represented in figure D-2. For notation purposes, k; is the wave number in the re- 
gion 1 (i=1,2,3,4). In regions 2 and 3, the s and d subscripts indicated the wave 
number is for the shear (transverse) and longitudinal (dilatation) waves respective- 
ly. 

In region 1 (exterior fluid), the total pressure, P1(r,@), is the sum of the 


incident plane wave and scattered pressure 
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Pi(r,8) = Po da En cos(nd) ( by Hp (kyr) + Jp(kyr)). (D-26) 


n=0 


In the viscoelastic layer, region 2, the potentials are 


ga(r,@) = Po x En cos(n@) (Cydn(kgar) + dnYn(kg2r)) (D-27) 
n=0 

W2i7,?) = Po) sn En Sin(n@) (endJn(ksar) + fnY¥n(ks2r)) . (D-28) 
n=0 


In the elastic layer, region 3, the potential functions are 


$3(r,8) = Po x En COS(N8) (GnJn(kg3r) + AnY¥n(ka3r)) (D-29) 
n=0 
W3(r,@) = Po) at En Sin(n@) (Indn(Ks3r) + mpYn(ks3r)). (D-30) 


n=0 





In the interior fluid, region 4, the pressure is written as 


© 


Pa(r,@) = Po) u" €, cos(n@) ( Ny Jp(kar)). 


n=0 


The displacements in terms of the potential functions are 


0 ] oF 
Te EN =e r 06 
1 @ ay 
aR EUE? =, : ar 


The stresses in terms of the potential functions are 





2 2 
T= od K 6 + 2p 3 : : = he on 
or Yr og r 06 or 
2 a 2. a ae OC«<‘( Ce 
re.ti‘ 2 ey oe a 
Yr rel] Yr 06 or or r or 


| 
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(D-31) 


(D-32) 


(De33) 


(D-35) 


The radial displacement in regions 1] and 4 may be written in terms of 


the pressure using Euler’s equation 


1 a 





(D-36) 
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The radial stress is the negative of the pressure in these regions. 


The boundary conditions are the requirement of continuity of stress 


and displacement across the three interfaces. They are as follows: 


(i) at re=c, 


] ] 2 2 
er ee ee Foes 
(ii) at re=b, 
2 3 ] 2 2 ee 3 
Ta ee els ae tl VO SOV = oF 
rr rr ré ré 
(111) at r=a, 
3 3 4 3 
yee = ler ; ele T 9 =a (D-37) 


A set of ten equations in ten unknown coefficients results as follows: 


Dx =A (D-38) 
n n 
where x is a column vector of the unknown coefficients, ahs is a col- 


and A,, and D is a 10 x 10 


umn vector of the nonvanishing elements A, , 


matrix with 60 nonzero elements which results from application of the 


boundary conditions. The D matrix nonzero coefficients are: 





2 
Py 
Ho 
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(2n° - Ke of “yd, (k 


2 


oe) 


2k ocd (kK aoe) 


2 2 
(Crm calnikancle= 2k cv ik. °C) 


en(k cd (k.4c) 


= en(k cY (k.4Cc) 


=k.,c Y’ ,k 


I, 
-kyc HO (k)¢) 


K ACs: , 6k 


d2 42°) 


d2 42° 


=n J ikea) 


- en(¥ (k 


n ¥ (ko) 


2n(I (kyo) - 
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: J bkooc)) 


SS) 


Ee cUm(kyec))) 


Ky nae Kao) ) 


; Le 2 
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eke) 
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cY (kc) 


Ka go?) 


znd 
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=n y 
n 


“Ka3 


5 bbe 


= -ny 
n 

= -nd 
n 

-ny 
n 
sg Soe 


Sb ee) 


nJ 


- K 9b Yk yo) 


(K.5b) 
(k.) 
bJ (kb) 


bY (k 445) 


= - nd (kb) 


(kK 35) 
(k 45D) 


(ka) 


asbeEe 
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+ (kc 


zg 
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- 2n Ke) 
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= ny (kb) 


ae sea 
ROT Reg) 
2 aes ((2n Ke b*)a (Kee yeceek. bd (ki b)) 
HM $2 n- d2 d2 n- d2 
= ce ((2n° 2 K2 b°)Y (k,.b) - 2k,,bY’ 6k b)) 
M. SZ n‘ d2 d2 d2 
"2 
= “hy en(k bd (kb) 2 J (k.,b)) 
a, 
. Ty an(k_ bY’ (kb) : ¥ (kb) 
| 2 ‘ 2 
= 2k ,bJ'(k4.b) - (2m - kod )d (kb) 
2 ‘ 2 
= 2k bY’ (k,.b) - (2n” - k 2b )Y (kb) 


én(J_(k.,b) - k.,bd’(k_.b)) 


2n(Y¥ (k 4b) - k 4bY’(k .b)) 


2n(J (kod) - kyobd’ (kb) ) 





a 
tt 


© 
1 


— 
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| 


© 
tt 


© 
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© 
tt 


2n(Y (kKyod) - kyobY’ (k4ob)) 


$2 


2 
= 2k bY’ (k 4b) + (ko ab 


See De DU (ike,_b) 


oe 


Bo 


B3 


ae (2k .bJ (kb) 


2 


3 
5 (ee OS 3) 


(2n fe. Kena a )J Kaa) 


(2n_ - Ke ya )Y, (k 


43?) 


AL) J 


2 2 
= 2k ,b0’(k_ 4b) + (keyb 
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2 
- 2n°)d_ (kb) 
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- 2n°)¥ (kb) 


; Jyh gg0)) 


—— an(k bY’ (kb) - 


=e 
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Tredse! 


(ke, ae 2n°)J(k_,b)) 


2 2 
(Koby - 2n")Y (k_4b)) 


- 2k aJ (ka) 


- 2k,.aY’(k 
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d3 


d3 43?) 
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vA 


fa gt eee oe 


oO 
t 
= 
Cc, 
= 
~-_ 
eam 


D =n ¥ bko3a) 

D = “kyad) (ka) 

D10.6 : 2n(J_ (kia) - k gad’ (ka) ) 
DI 0,7 - en(Y (ka) - K ga’ (Ka) 


2 


; 2 2 
D = 2k ad (k 4a) + (k 44 - 2n )J (ka) 


2 


2 2 
D = 2k ,aY’(k.4a) + (koga - 2n°)Y (k 3a) 


The A vector has two nonzero elements which are 


2 Kc 


c 
A, = - a J Ck, ¢) and A. = ° 
2 pie 





J’ (ky¢) (D-40) 


The D matrix has been adjusted by multiplying column 1 by pyu and column 10 by 
paw as suggested by Gaunaurd>?, (The above coefficients were originally derived 
by Guanaurd>? but were rederived for verification. Verification of the coefficients 


derived by Flax and Neubauer>® could not be obtained.) Additional scaling has 
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been done on the sixth and seventh columns to prevent computational problems. 
For the n=0 case, the ¥2 and ¥3 potential functions are 0 due to the sine 
term. This results in 7;g always being zero regardless of the constants, which 
causes these three boundary conditions to be meaningless. In addition, the conti- 
nuity of v displacement is also meaningless. The overall result is to reduce the D 
matrix from a 10 x 10 matrix to a 6 x 6 matrix by eliminating rows 3,5,7 and 10 and 


columns 4,5,8 and 9. 


Elastic Shell Derivation 


The solution for the scattered pressure from the purely elastic shell is 
obtained from the preceding discussion by eliminating the viscoelastic layer (layer 
2) and applying the outer fluid boundary conditions to the elastic layer. The total 
pressure in fluid 1 is expressed by equation (D-26). The vector and scalar 
potentials for the elastic solid are expressed by equations (D-29) and (D-30). The 
total pressure in the interior fluid (fluid 4) is expressed by equation (D-31). The 
boundary conditions are: 


(i) at r=b, 


(iii) at re=a, 





We 


= -P U =u, Tf , = 0 (D-41) 

A set of six equations in six unknown coefficients results as follows: 
D Naps A (D-42) 

where x is a column vector of the unknown coefficients, A. is a col- 


and A., and D is a 6 x 6 


umn vector of the nonvanishing elements A, 2 


matrix with 28 nonzero elements which results from application of the 


boundary conditions. The D matrix nonzero coefficients are: 
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~én(k j.aY (kia) - ¥ bka34)) 


5,3 
ea ea tana- eanayor (koa) 
5,4 s3 on $3 $3 noass 
Deg = 2k. ,a¥/ (koa) + (Ke 4a - ane )Y_(k.9a)) 
De 2 = K 3 a J’ a Ky3a) 
eno doen ni das 
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The A vector has two nonzero elements which are 


v kb 
A= -—— 9) (kyb) and A, = 3 JN(k,b) 


3 pw 


(9S 


(D-44) 


As with the layer solution, the D matrix is reduced for the n=0 case (D matrix 


reduces to a 4x4 matrix). 


The elastic shell code is verified by comparing the theoretical results with the 


experimentally obtained form function for the simulated infinite cylindrical shell. 


In addition, the theoretical computer results were qualitatively compared with the 


theoretical form functions illustrated by Flax and Neubauer-®. The viscoelastic 





176 


layer program was verified by: (1) recovering the experimental shell form function 
when the viscoelastic layer parameters simulated a lossless water layer, and (2) the 
pressure release form function generated in Appendix A was recovered when the 
viscoelastic parameters were Set to simulate a low loss air layer. 

The Bessel functions are generated by the methods indicated in Appendix A. 
The accuracy of the Bessel functions were checked against published tables and 
other library routines. In addition, the Wronskian values calculated using the 
Bessel functions and their derivatives were continually compared against the exact 
values for a estimate of solution accuracy. The matrix equations are solved by a 
standard IMSL routine (LEQ2C) which utilizes a lower-upper decomposition 
method with iterative improvement. All computations were performed using 


double precision on a Vax-750 or Micro-Vax computers. 
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APPENDIX E 


The Treatment of Multiple Reflections in the Panel for 
Plane Wave Modulus Experiment 


Due to the large difference in the specific acoustic impedances of the water 
and the materials used in the plane wave modulus experiment, significant 
reflection will occur at all interfaces. The purpose of the experiment is to 
determine the transmission of sound through the material. The reflections of im- 
portance are the internal reflections that occur in the material which eventually 
add to the directly transmitted pulse but with delayed phase and reduced ampli- 
tude. To fully assess the transmission through the panel, these multiple reflections 
must be taken into account. For the closed cell rubber materials, the attenuation 
in the material is so large that the multiple reflections are significantly reduced and 
hence have little or no effect on the direct transmission. The corprene and cork-ni- 
trile composites, however, do not have a large enough attenuation to permit 
neglect of multiple reflections. 

As illustrated is figure E-1, an incident wave contacts the surface (2) and a 
portion of the wave is transmitted. At the next interface between the two materials 
a portion of the first transmitted wave is transmitted into the medium. The portion 


transmitted from the water into the material is related by the transmission 
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Figure E-1 Multiple Internal Reflection in a One-Dimensional Panel 
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coefficient Tj 2 and the portion transmitted from the material into the water is re- 


lated to Tz 4 where 


and, rj is the specific acoustic impedance of the water and ro is the specific 
acouStic impedance of the panel material. There is also reflection of the sound 
wave occurring at the material-water interface with the amount of reflection 


determined by the reflection coefficient, R, where 


= (E-2) 


The incident wave, Po, is related to the transmitted wave, P;, by 


_ -al twl/c 2 -3al t3w/c 
P. = Py { Maloy e oe + Ty ato yh e e 


| earns eee! (E-3) 


2 oat at wl/c P 
1 2e2.1 ae 


The direct transmission can be factored out resulting in 
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po 2 a4l _rwal/c oat } (E-4) 
een MS tect a deen ie) - CES then 
Meee 1 
equation (E-3) can be written as 
p 
SCNT s Wax): (E-5) 
0 
Since x<]l, the series sums to = Equation (E-5) becomes 
P 
t T 
am S oes a 


From the experimental setup, the ratio of transmitted pressure 


to incident pressure is 


wo 


t ae PA) | bol 
ae x0; |= 


ie exp( Ld ) oe 


~~ 


where ee is the unwrapped phase for the water only, bo is the un- 


wrapped phase with the panel in place, C is the propagation speed of 
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the water and c is the phase speed of the sample. The ratio of amp- 


litudes, A, and phase, ¢, is determined directly from the FFT data with 





Be” Earp OO set Sa (E=/) 


The — term is due to a thickness of water, L, being displaced by 
vs 
the panel. 


Iteration Process 


The equation to be iterated is 
T = A exp(tg) (1-x) (E-8) 


1. No internal reflections is assumed, i.e. x = 0. 





the initial phase speed is c = (E-9) 


From the initial phase speed, 1 0 and 1, , are calculated 


and the initial attenuation, a, is determined from 


A= No M4 exp(-aL). (E-10) 


2. From the initial values of c and a, the phase and amplitude 


corrections for the multiple internal reflections are cal- 


culated. 
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